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PREFACE 



This is one of a continuing series of reports of the Ford Foundation 
sponsored Research Program in University Administration at the University 
of California, Berkeley. The guiding purpose of the Program is to under- 
take quantitative research which will assist university administrators 
and other individuals seriously concerned with the management of university 
systems both to understand the basic functions of their complex systems 
and to utilize effectively the tools of modern management in the alloca^ 
tion of educational resources. 

In this paper we present SPACE> a Space Planning And Cost Estimating 
simulation model designed to allow analysis of alternative class scheduling 
patterns and their consequent resource demands. Several illustrative 
examples of the model^s use. are given, with documentation of the valida- 
tion procedures using data from the University of California. The examples 
show that changes in class scheduling patterns have no significant effect 
on total operating and capital costs. A listing of the computer program 
and the input data specifications are included in the Appendix. 

This research was initiated and conducted in the Office of the Vice 
President — Planning of the University of California and was partially 
supported by the California Coordinating Council for Higher Education. 
The authors wish to express their gratitude to Leonard Jay Jacobson of 
Mathematica for patient counsel which resulted in the avoidance of many 
pitfalls; to Maggie Hardy; and to John Lafler, whose many insights and 
ingenious "mathemachinations" gifeatly improved the initial design, and 
whose programming of the FAMSIX module made the whole thing possible. 
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INTRODUCTION 

A 1970 resolution by the California legislature called for an extra- 
ordinary increase in classroom utilization by the University of California 
and the California State Colleges. The California Coordinating Council 
for Higher Education (CCHE) sought to develop an objective mechanism for 
estab?^ishing realistic utilization standards with appropriate consideration 
of any demonstrable effects on both operating and capital costs. Mathematica, 
Inc. , undertook and completed a iiiammoth mathematical modelling effort under 
severe time constraints."^ Unfortunately, the Facilities Analysis Model 
(FAM) was not sufficiently flexible to allow detailed analysis of the 
alternatives under consideration. The SPACE project began as an attempt to 
modify FAM to be more accurate and more useful to physical planners in 
higher education. 

The University's efforts at modification were soon abandoned as un- 
workable, and a new model called SPACE (Space Planning And Cost Estimating) 
was developed, drawing heavily on Mathematica 's experiences but differing 
substantially in several fundamental respects. As result the SPACE model 
is in operation at the University and, with a few modifications, is being 
prepared for implementation by the California State University and Colleges' 
Office of Analytic Studies. 

SPACE is a computerized simulation model designed to provide planning 
data for managers of institutions of higher education. It is primarily 
concerned with the physical facilities for departments of instruction and 
research. SPACE relates enrollments, clasc scheduling policies, class 

Hlathematica, Inc., ''The CCHE Facilities Analysis Model," Vols. I-IV, 
Princeton, New Jersey, December 1970. 
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size policies, and faculty workload policies to physical resource needs 
and the resulting utilization rates of classrooms and class laboratories. 
SPACE also computes the faculty needs, the associated current costs, and 
converts the capital costs to comparable annual costs. 

In its present stage of development, SPACE is confined to the '*I&R" 
departments and thus cannot serve all the planning needs met by RRPM or 
CAMPUS. On the other hand, its moderate size and its cone jition on 
physical resources provide an easily manageable tool for answ^.: Jng a sig- 
nificant number of what -if questions at a level of detail not available 
in the more general models. Furthermore, to our knowledge, SPACE is the 
only operational model with the capability of simulating both the time 
scheduling of classes and the assignment of classes to instructional rooms* 
As a tool for University planning SPACE can help improve the physical 
planning process by augmenting and refining the calculations uf facilities 
needs — calculations frequently mad^ without e:j:plicit consideration of such 
factors as class size policies and faculty workload policies* 

It is with regard to these functional policy alternatives that SPACE 
should prove most useful. Typical what- if questions within the domain of 
SPACE include: 

1. What if class schedules are extended beyond the normal daytime 
hours? 

2. What if faculty workloads and/or staffing patterns are changed? 

3. What if a specific new program is added (or a specific existing 
program is eliminated) at a particular campus? 

4. What if classes larger (and/or smaller) than certain sizes are 
eliminated in certain (or all) programs? 

Answering such questions in terms of the probable effects on resource 
requirements is a fundamental role of both statewide and campus academic 
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and facilities planners, SPACE can reduce the staff effort and time 
required to perform these calculations and should also substantially 
improve the accuracy of the answers* 



> 



MODEL DESCRIPTION 



Model Overview 

SPACE is composed of two separate modules, FAMSIX and VARICOS, each 
of which is written in FORTRAN IV. The modules can be run separately or 
in tandem, with FAMSIX providing a machine-readable file for use as input 
to VARICOS. 

FAMSIX performs sequentially the 6 basic operations lasted below. 
These are repeated for each yiar of the planning horizon. In its current 
form, the maximum number of years in the planning horizon is 10. 

1. A Gaussian pseudo-random number 2, -aerator operates on /eekly 
student hour input to produce preliminary total of weekly 
room hours* The level of aggregation is class size range, 
discipline, class type (lab or non-lab), and course level. 

2. Classes are assigned to hours of the i^eok according to a 
priority schedule which reflects the user*s preferences for the 
distribution of class hours. (Note: Classes are treated as 
integral units in that no fractional class hours are generated.) 

3. The effects of scheduling on average class size (time-preferences) 
are computed using a combination of linear and exponential rela- 
tionships. The initial distribution of classes to class size 
ranges is adjusted to account for this phenomenon while keeping 
the original total of weekly student hours intact. 

4. Classes of each size are assigned to rooms of appropriate sizes 
by a randomization routine governed by user-defined parameters. 

5. For each hour of the week, the current inventory of classrooms 
and class labs is compared to needs. If rooms are needed but 
not available, the inventory is checked for the availability of 
rooms which are larger than required and the overflow is assigned 
where possible. If there are no larger rooms available, or if 
there is an inadequate number, rooms of the appropriate type and 
size are "constructed.'* 
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This is one of a continuing series of reports of the Ford Foundation 
sponsored Research Program in University Administration at the University 
of California, Berkeley. The guiding purpose of the Program is to under- 
take quantitative research which will assist university administrators 
and other individuals seriously concerned with the management of university 
systems both to understand the basic iunctions of their complex systems 
and to uf'^.lize effectively the tools of modern n^anagement in the alloca- 
tion of educational resources. 

In this paper we present SPACE, a Space Planning And Cost Estimating 
simulation model designed to allow analysis of alternative class scheduling 
patterns and their consequent resource demands. Several illustrative 
examples of the model *s use are given, with documentation of the valida- 
tion procedures using data from the University of California. The examples 
show that changes in class scheduling patterns have no significant effect 
on total operating and capital costs. A listing of the computer program 
and the input data specifications are included in the Appendix. 

This research was initiated and conducted in the Office of the Vice 
President — Planning of the University of California and was partially 
supported by the California Coordinating Council for Higher Education. 
The authors wish to express their gratitude to Leonard Jay Jacobson of 
Mathematica for patient counsel which resulted in the avoidance of many 
pitfalls; to Maggie Hardy; and to John Lafler, whose many insights and 
ingenious "ma thema china tione" greatly improved the initial design, and 
whose programming of the FAMSIX module made the whole thing possible* 
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6. Utilization statistics are compured for each classroom size and 

for each class lab size, the class labs being further distinguished 
by discipline and level. 

VARICOS is basically an accounting module with 4 component operations 
which are repeated for each year of the planning horizon: 

1. Weekly room hours (generally but not necessarily from FAMSIX) and 
enrollment projections are operated on by functions defined in 
terms of workload and other parameters to compute faculty staff- 
ing needs by rank. 

2. Faculty office and research space needs are computed along with 
the associated support space needs. These needs are compared to 
the existing inventory and new space is "built" where needed, 

3. Faculty salaries, support costs, and facilities maintenance and 
operating costs are computed. 

4. Two kinds of capital costs are computed. The first are project 
costs, i.e., what is commonly called total capital outlay. The 
second set of costs are the annual, cash flow needed from the 
State or other funding agency. It is assumed that physical faci- 
lities are financed by the Issuance of bonds, and that the issues 
are sequenced according to the costs ii.curred during the planning 
and construction period (i.e., the 4 years prior to the first 
usage of the facility). 



FAMSIX - Module Description - 

This section details the operating logic of the six components of the 
FAMSIX module described briefly above. These calculations are repeated 
for each iteration of the test run. A program flow chart is included to 
assist the reader in understanding the sequencing of the components and 
their interrelationships within the model. The underlying mathematics are 
Q not spelled out, but examples are provided to give the reader a feel for 

ERIC 
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the processes involved. An attempt has beer, made to provide enough de- 
tail so that a progrannner examining the program code will have adequate 
guidelines but not so much that the non-technician will be bogged 'lown. 

The Generation of Initial Weekly Class Hour Totals from Weekly Student 
Hour Inputs 

FAMSIX is driven by Lab and Non-lab Weekly Student Hour inputs. 
(These are displayed at the user's option in Table 101.) Historically 
the distribution of WSH for any given discipline and course level varies 
greatly across class sizes, both for Non-lab and for Lab classes. Since 
it is Impossible to predict with certainty these distributions, FAMSIX 
operates under the assumption that the proportion of WSH in each of the 
15 class size ranges varies normally for any category. It further assumes 
that the average class size in each class size range is constant. 

Thus, for each discipline, course level, and class type, FAMSIX cal-- 
culates the proportion of WSH in each class size range by a Gaussian pseudo- 
random number generator. For each category, the random number generator 
is given the mean proportion and the standard deviation tor each class 
size range. It returns the proportion of WSH in the class size range (set 
equal to zero if negative). Since the 15 proportions thus determined 
generally will not sum to 1.0, they are scaled proportionately until they 
do sum to 1.0. The resulting proportions are then multiplied by the total 
WSH for the category to determine the number of WSH in each class size 
range. 

Each of these numbers is then divided by the (fixed) average class 
size in the respective class size range to give initial numbers of Weekly 
Class Hours. Since the average class size generally will not exactly 



7 



divide the number of WSH, the resulting WCH figures are rounded. These 
"integerized** WCH numbers are the numbers on which FAMSIX continues to 
operate. They are not the final numbers of Weekly Room Hours. They are 
displayed at the user's option in Table 211 (Non-Lab) and in Table 212 
(Lab). 

Note that if desired the mean values for class size distribution 
may be chosen to be negative inputs. In such cases, the category will 
have extremely different class size distributions from one iteration 
(year) to the next, and may even disappear in some iterations. This is 
a useful device for a category (discipline, level, and class type) in 
which such extreme variation is characteristic, either in fact or because 
the data are reported inconsistently. 

The Assignment of Weekly Class Hours to Hours of the Week 

The assignment of Weekly Class Hours to the hours of the week is made 
through the mechanism of a "priority schedule." There is a separate 
schedule for each class type (Lab and Non-lab). The priority schedule 
associates with each hour of the week zero, one, or more priority numbers. 
The relative number of such priority numbers associated with each hour 
represents the relative distribution of class hours to the hours of the 
week. The size of the priority numbers associated with any hour reflects 
the degree of preference for classes in that hour. Thus, if the 37th hour 
is assigned one priority number and the Ist hour is assigned 5 priority 
numbers, class hours will be scheduled into the 1st hour 5 times as often 
as they are scheduled into the 37th hour. Furthermore, if the priority 
number assigned to the 37th hour is large, e.g., 110, and if the set of 
priority numbers assigned to the 1st hour includes smaller values, 
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e.g., 1, 15, 70, 105, and 116, then the 1st hour will tend to have larger 
classes scheduled in it than WlII the S'^th hour. 

The medhanism works like this: There are (in the model) 96 possible 
hours of operation in a week. There are 250 possible priority numbers, 
p , of which n are chosen to be associated with the 96 hours (p = 1, 2,.. 
n <^ 250). FAMSIX treats the class hours as if they were ordered in a long 
list beginning with the largest classes and ending with the smallest* With 
i.n each class size range the class hours are ordered by course level and 
within course level by discipline. Then for hour h , with priority n'omber 
p , the p^^ class hour in the list and every n^^ class hour in the list 
thereafter are scheduled for hour h . In the example above, if n were 
200, hour 1 would receive the 1st, 201st, 401st,... class hours in the list 
It would also receive the 15th, 215th, 415th,...; the 70th, 270th, 470th,.. 
the 105th, 305th, 505th,...; and finally the 116th, 316th, 5l6th,.*. class 
hours in the list. Hour 37 on the other hand would receive only the 110th, 
310th, 510th, class hours. Clearly the 37th hour will receive approxi- 
mately one-fifth as many class hours as the 1st hour. Because of the 
ordering of the list froir large classes to small, the 37th hour will also 
receive smaller classes on average than the 1st. 

The Effects of Scheduling on Average Class Size - "Hour-Size Function" 

Historical data gathered at the University reveal that, ceteris 
pcLTibuSi class size varies according to the hour at which the class is 
scheduled and the number of classes which are scheduled in the hour. Most 
classes would attract fewer students at 6 p.m. than at 10 a.m., and fewer 
on Saturday than on Monday. Furthermora, it is obvious that if a great 
many classes are offered in any hour, the last one scheduled will be 
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smalleL than it would have ha-.in if it had been scheduled at an equally 
attractive hour in which fewer classes were scheduled. This is true 
because some of the students who would otherwise be attracted to the class 
are more likely to have enrollea in other classes offered at that hour. 

FAMSIX simulates this effect by altering the average class size in 
any given hour according to the *'hour-size function," This function relates 
the expected average class size to an altered average class size by multi- 
plying the expected size by a factor f , The product is the altered 
average class size. The factor f is determined by the following function: 

f = c • p^; 

where: 

c = a scalar multiplier which is a function of the hour in question, 

p = the proportion of WSH in the given hour, and 

y = some exponent,^ 
The application of the function to the Weekly Class Hours generated by 
FAMSIX is somewhat complicated since FAMSIX assumes that the average class 
size in each class size range is fixed. Furthermore, it is not desirable 
to increase or decrease the total number of WSH appreciably. Consequently, 



The problem of determining reasonable values for c and y can be 
greatly simplified Ly the following procedure: 

1) Set y equal to some value near -0.25* This value has been 
shown to be reasonable for U,C. studies. 

2) Determine the approximate proportion of weekly student hours 
in each hour of the week according to the way classes are 
actually scheduled, 

3) Solve 

f = 1*0 » c p^ 

for c for each hour. For the actual schedule, the function 
should leave class size unchanged. (Presumably the class-size- 
distribution inputs are derived from actual dataO Thus, using 
actual proportions, cp^ should equal 1,0, 

O 

ERLC 
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FAMSIX leaves the number of WSli in each hour intact and changes the distri- 
bution of classes to the class size ranges. The net effect is that the 
number of class hours is altered so that the average class size over all 
class size ranges is increased or decreased by the factor f . This new 
number of class hours is then distributed to the disciplines and levels in 
the same proportions as the original number of class hours. Care is taken 
so that the number of weekly student hours within a discipline and level 
is not altered appreciably by this process. 



Exampl e 

Suppose that there are only 3 class size ranges and that the fixed 
average sizes in those ranges are 10, 15, and 20. Further assume 
that for the hour in question the factor f = 1.15, and the initial 
distribution of class hours is as in the table below. 

Size 1 Size Z Size 3 Weekly Student Average Class 
(10) (15) (20) Hours Size 

15 20 0 450 12.86 



FAMSIX first changes the 20 classes of size 2 (i.e., 15 students each) 
so that the average will be, as nearly a3 possible, 15 x 1.15 = 1/.25, 
re*^.aining closely the 300 WSH generated by the original 20 classes of 
15 students each. 

Size 1 Size 2 Size 3 WSH ACS 

9 8 295 17.35 



FAMSIX then notes that it is 5 WSH short. Adding these to the 150 
generated by the size 1 classes, it changes these 15 class hours so 
that the average will be nearly 10 x 1.15 = 11.5. 

Size 1 Size 2 Size 3 WSH ACS 

9 4 ,150 11.5 



Nothing could be done with the 5 WSH, so they are lost. (The greatest 
error in this process is a gain or loss of 1/2 the fixed average class 
size of the smallest class size range.) 

The net effect of these changes is summarized below. 
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Size 1 Size 2 Size 3 WSH ACS 

9 13 8 445 14.8 



The altered average class size target was 1J5 x 12.86 = 14.8 • 
Note that the initial number of class hours is 35. The final number 
determined hy the hour-size function is 450/14.8 - 30 . 

Assignment of Classes to Rooms 

Each quarter or semester, some time before students enroll in their 
courses, departments determine their course offerings and the campus 
Registrar (or other class scheduling agency) must determine the facilities 
to be reserved for those courses. Using departmental predictions and 
past experience as guides, the Registrar is able to make some estimates of 
the class sizes which will be generated by those courses. Since it is 
generally impossible to determine the exact course enrollments and class 
sizes, the distribution of actual classes to class size ranges will differ 
from the predicted distribution. In an effort to minimize the amount of 
facilities reassignment, tha Registrar will generally try to provide rooms 
so that classes which turn out to be larger than expected will seldom have 
to be reassigned to other facilities. The effect of this hedge against 
possible chaos is some "looseness" in the relationship of the sizes of 
the classes to the sizes of the rooms to which those classes are assigned. 

FAMSIX simulates this effect in two stages. The first stage is accom- 
plished with the help of the pseudo--random number generator. For each 
class size range, the probability distribution of classes to room size 
ranges is an input to the model. Each class hour generated via the hour- 
size, function routine is associated with a "request" for a room (classroom 
in the case of non-lab classes and class lab in the case of lab classes.) 
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The size of the room is determined by the randoDi number generator operating 
on the distribution parameters- 

Each class hour is treated separately so that fractions of classes 
are not associated with toom ^'requests." For example, suppose: 

0.25 * probability that class of size 3 is associated 
with room of size 3; 

0-75 probability that class of size 3 is associated 
with room of size 4; and 

5 = number of class hours of size 3 generated in 
some hour. 

A strict application of the probabilities world yield 1,25 classes assigned 
to rooms of size 3, and 3.75 classes assigned to rooms of size 4, Since 
FAHSIX treats class hours as integral units, such a method is inadequate. 
FAMSIX, in fact, could assign all 5 classes to rooms of size 4, ^r all 5 
to rooms of size 3- On average, however, 1/4 of all classes (across all 
hours and years) of size 3 would be associated with requests for rooms of 
size 3 and the rest with rooms of size 4- 

For reasons which will become obvious after reading the next section, 
the parameters of the class-size-to-room-size distributions should not be 
identical to historical distributions. Rather, they should reflect the 
Registrar's effort to avoid chaos in logically identifying room size re- 
quests with predictions of class sizes. Thus, if it were possible to know 
class sizes with certainty before making room assignments, the most reasonable 
distribution would assign a class of size x to a room of size x with 
probability 1. 



13 



Spill -Up and New Construction 

The second stage of the class to room assignment routine is the 
"spill'-up" operation. For each hour of the week the room requests are 
compared with the current inventory, beginning with the requests for the 
largest rooms. If the Inventory is inadequate to the needs for a particu- 
lar room size, rooms of larger sizes are made available if they have not 
already been claimed. If rooms are still needed after the search has been 
made for unclaimed rooms of that size or larger, the model "builds" them, 
and they are added to the inventory. It is this spill-up process that 
requires the class-size-to-room-size distribution discussed above to be 
a theoretical distribution. Some spilling up occurs in fact and would be 
reflected in a historically derived distribution. Using the historical 
distribution would, therefore, result in spilling^up twice. 

FAMSIX has one category of classrooms, and thus assumes that any 
non-lab class can be assigned to any classroom, provided only that the 
room is large enough. This assumption is not, of course, entirely realiS'- 
tic because some classrooms should have special facilities for certain 
non-lab classes — e.g., lecture-demonstration facilities for certain science 
lecture classes — and those specific classes should be assigned to those 
rooms, even if some of those classes are much smaller than the room. But 
it is Impractical to include the capability of simulating those and other 
realistic complexities in FAMSIX*s handling of non^lab classes, not only 
because the computer core requirements would be greatly increased, but 
also because the input data requirements would be inordinately large. 
Presumably, however > this simplification does not cause the classroom 
requirements computed by FAMSIX to be significantly different from those 
which would be generated by the more "realistic" model, at least in terms 
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of the total required floor area. 

Lab classes are treated differently than non--lab classes. A lab 
class must be assigned to a class lab appropriate to the particular level 
and discipline of the class. Thus, all Arts lab classes arc assi^.gned to 
Arts class labs. Currently, there is a similar one-to-one relationship 
for the level of the class and the level of the laboratory. This identity 
is a proxy for the real-life non-interchangeability of labs between sub- 
disciplines of a given discipline and some non-interchangeability among 
levels. Future versions of the model will allow lab classes of level x 
to be assigned to class labs of level y , providing that level y is at 
least as high as level x . 

The final distribution of classes by size to rooms by size for each 
hour is displayed in Tables 401 (Non-lab) and 402 (Lab) at the user's 
option. 

Utilization Statistics and VARICOS Interface 

After all preceding steps are completed for all hours of the week, 
summary data are computed. Table 121 displays the final totals of WSH 
by level, discipline, and class type. These are computed by multiplying 
the final numbers of WCH (now identified as Weekly Room House, WRH) in 
each class size range by the fixed average class size in that range. 

The final numbers of WRH for each discipline and level by class size 
are displayed in Tables 221 (Non-Lab) and 222 (Lab). The summary of these 
data, excluding class size information, is stored externally for use by 
VARICOS. The construction necessary to meet room requests is displayed 
in Tables 501 (Non-lab) and 502 (Lab). The classrooms to be constructed 
for use in the year in question are distinguished by size only. The labs 
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WCH by size VIA (I) 

HOUR-SIZE FUNCTION 



ASSIGN CLASSES TO 
ROOM SIZES 



(J) 



COMPARE ROOM NEEb3 
TO INVENTORY; SPILL 
UP WHERE POSSIBLE; 
BUILD WHERE NECESSARY 




COMPUTE NEW WSH 
TOTALS & SUMMARIZE (L) 



FINAL WRH 



OTALS 




\221 & 222 



UPDATE INVENTORY 
ASF 



COMPUTE UTILIZATION.MN 
STATISTICS 




JABLES 301 
& 302 



TABLES 501 
& 502 I — ^ 



WRH BY 
DISC, LVL 
FOR VARIC OjJ 



NEW ASF & 
INVENTORY 
BY DISC FOR 
VARICOS 



♦Letters (A) - (N) are referenced in program listing. See Appendix. 



16 



are further distinguished by discipline and level. 

The inventory of assignable square feet is calculated along with the 
ASF of newly completed space. This information is tagged by level and 
discipline for class labs. Classrooms are arbifiarily identified with 
Discipline #1. The ASF data c^re also stored externally for use by VARICOS. 

Tables 301 (Non-lab) and 302 (Lab) are the final outputs of FAMSIX. 
They are displays of sunmiary utilization rates for each room size and in 
total. (The information is further disaggregated into disciplines and 
levels for laboratory utilization summaries). The information displayed 
separately for classrooms and class labs includes the current inventory 
of rooms and stations, the number of WRH, WRH per room, Unweighted Station 
Occupancy, WSH per station, ASF per station, and ASF per WSH. 

VARICOS " Module Description 

This section details the operating logic of the four components of 
the VARICOS module described b-'iefly in the model overview. These compo- 
nents are repeated for each iteration (year) of the tesc run. The general 
level of detail corresponds to that of the FAMSIX module description. 

Faculty Requirements 

Faculty needs are generated from enrollment projections. Lab and Non- 
lab weekly class hours (i.e., the weekly room^-hours determined by FAMSIX), 
and a combination of stijdent'-f acuity ratios and wrkload parameters. The 
iQutine provides for special handling of Teaching Assistants because of 
the Important instructional role played by T.A.'s at the University of 
California. 
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VARICOS computes? f'^vr categories of faculty: T,A.'s and 3 categories 
(not necessarily ranks) of rion-T.A. faculty. The non-T.A. faculty will be 
referred to as "regular** faculty in rhis description. The computation 
takes place in four steps. 

Step 1: Dividing the Vorkload 

FAMSIX provides VARICOS with the number of WRH in each class type, 
course level, and discipline. VARICOS inputs include the proportion 
of WRH taught by T.A.'s and the proportion taught by regular faculty 
in each category of WRH. These proportions are applied to the appro- 
priate WRH figures to give the number of WRH in each category which 
are taught by T.A.'s and the number which are taught by the regular 
faculty. 

Step 2: Assigning FTE to the Workload 

VARICOS inputs also include the number of Weekly Faculty Contact 
Hours (WFCH) p'^r FTE, separately for T.A.'s and regular faculty, in 
each of the WRH categories. The division of T. A. -taught WRH by the 
workload parameter (WFCH/FTE-TA) in each category gives the number 
of T.A."s needed. The same operation using regular faculty workload 
parameters provides a partial tota?. of regular faculty needed. 
Additional regular faculty can be generated on the basis of T.A. 
supervision requirements. For each discipline and class type the 
supervision factor, regular faculty FTE per T.A. FTE» is multiplied 
by the number of T.A."s generated and the product is then added to 
the regular faculty subtotal. 

Step 3: Generating Faculty via Enrollments 



Additional regular faculty may be generated from student-faculty 
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ratios for each of the four student levels. Enrollments by level and 
major are multiplied by these ratios and the products added co the 
subtotals from Step 2. Step 3 thus serves two functions: first, it 
may be used as the basis for generating faculty strictly from enroll- 
ments by setting the workload parameters to zero; and, secondly, it 
may be used to augment the faculty generated from workload factors 
on the assumption that students generate a need for faculty to super- 
vise independent study or other work not reflected in the WRH figures. 

Step 4; Distributing the Regular Faculty to Ranks 

For each discipline, the proportions of regular faculty in each of 
the 3 ranks (or other categories) are data inputs. These proportions 
are applied to the regular faculty discipline total to give, with the 
T.A. figures, the final numbers of faculty in each rank and discipline. 

Office and Research Space 

There are 5 types of space requirements computed by VARICOS and FAMSIX: 
1) Classrooms; 2) Class labs; 3) Research labs; 4) Offices; and 5) Support 
space. The assignable square feet of the inventory and of the newly com- 
pleted space are compiled separately. FAMSIX provides the figures for the 
first two space types and VARICOS computes the figures for the last three. 
In each non~class type of space, space-standards parameters are used to 
determine the space needs. The needs are compared to the current inventory. 
If the inventory is inadequate, the additional space needed is noted and 
added to the inventory. Each discipline has its own space and hence its 
own standards. 

Research space requirements are computed as: 
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(ASF/Regular faculty) x (Regular faculty) + (ASF/TA) x (TA's) + 
(ASF/Level 3 majors) x (Level 3 majors) + (ASF/Level 4 majors) x 
(Level 4 majors). 

Normally, student majors of levels 3 and 4 are identified with first and 
second stage graduate students. 

Office space needs are computed according to the same formula, sub- 
stituting the appropriate space-standards parameters. 
Support space requirements are computed as: 

(ASF/Total Research & Office ASF) x (Total Research & Office ASF) 
+ (ASF/Total class lab ASF) x (Total class lab ASF). 

Faculty Salaries, Support, M&O-Plant 

Because FAMSIX and VARICOS are primarily intended to provide planning 
data related to cla^s scheduling, no attempt has been made to compute total 
operating costs • VARICOS therefore concentrates on those variable costs 
which are most directly related to scheduling and to class-size policies, 
viz., faculty salaries, departmental support, and physical plant maintenance 
and operation. 

Currently these computations are somewhat crude. Salaries are com- 
puted for each rank and discipline by multiplying the average salaries 
for each rank by the appropriate FTE figures. For each discipline the 
departmental support per regular faculty FTE is a data input. This param- 
eter is multiplied by the regular faculty FTE to get the support costs. 
There is a single M&O-Plant factor for all types of space. The total M&O 
cost is the product of this parameter and the total inventory ASF. 

All dollar figures are "constant salary dollars,'' i.e., faculty salaries 
and support costs are not inflated* M&O costs and construction costs are 
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inflated, however, since these costs hwe been rising more rapidly than 
faculty salaries. The difference between the two actual rates of infla- 
tion is the inflation factor applied to the space costs. 

Capital Costs and Outlay 

For each space type and discipline, the capital outlay per additional 
ASF is a data input. For each year of the run, these parameters are 
applied to the newly completed ASF to give the capital outlay generated. 
Since planning and construction generally occur over a 4-year period 
prior to occupation, the capital outlay is spread back over the four 
years prior to the year the facility is first added to the inventory. 

The outlay amounts, however, are not directly comparable to the salary, 
support, and M&O cost figures since they actually reflect a multi-period 
:ash flow from the funding agency. VARICOS assumes that all of the facili- 
ties construction is financed through the issue of bonds and that the 
issues are sequenced according to the pattern of cost incidence during 
the planning and construction period. Therefore, the outlay figures are 
distributed over the four--year period preceding the first occupation of 
the facility and bonds are issued to meet those costs as they are incurred. 
The debt-servicing of the bonds represents the economic costs of capital 
construction which are most appropriately compared to the operating costs 
already computed. It is the debt-service payments, therefore, which are 
displayed for each year of the test run, along with the operating costs, 
in Table C. For each iteration. Table B displays the year's capital outlay 
and current inventories. Similarly, for each iteration. Table A displays 
the year's faculty requirements and their associated costs. 



21 



MODEL TESTS AND RESULTS 

Model Characteristics 

All initial testing of the model was performed with data from the 
Santa Barbara campus of the University of California. When it was im- 
possible to generate parameters from the available UCSB data, data from 

3 

other general campuses of the University were used. The base year for 
all model tests was the 1969-70 academic year. 

It was found that SPACE is robust with regard to most parameter 
changes and Monte Carlo perturbations when viewed from a macro level. 
Overall class sizes, faculty requirements, facilities requirements, and 
costs vary only slightly in response to modest parameter changes. On 
the other hand, class size distributions and resource requirements vary 
radically within particular disciplines and levels. This kind of varia-- 
tion was not only anticipated It consciously sought because historical 
data show substantially greats : variation from year to year within dis^ 
ciplines than for the campus as a whole. It is one of the major goals of 
SPACE to demonstrate to University planners the ranges of needs which 
must be considered in facilities planning. 

Throughout the period of the model's development, numerous test runs 
were made to determine the sensitivity of the results of parameter changes, 
and to validate particular routines by comparing the outputs with histor-- 
ical data. Additional test runs were made to test the flexibility of 
the model in anowering what'-if questions. The reaults were uniformly 
favorable. In the three sections which follow, seven of the test runs, 
typical of the lot, are discussed in relation to these purposes. 

3 

Berkeley, Davis, Irvine, Los Angeles, Riverside, San Diego, Santa 
Barbara, and Santa Cruz. 
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Validation Tests; Comparisons with UCSB History and with Current Planning 
Methodology 

For the purpose cf validating the model, two types of tests were made: 

(1) straightforward comparisons of SPACE outputs with historical UCSB 

data, and (2) comparisons of SPACE-^enerated assignable square feet (ASF) 

requirements with the corresponding requirements computed by University of 

4 

California planning officers using the traditional methodology. 

The first type of validation test was made using historical UCSB data 
for the majority of SPACE inputs • Where UCSB data were incomplete or 
inadequate, data from the University's other general campuses were adapted. 
Enrollment data were estimated by assigning target three^-term-average 
headcount enrollment totals for undergraduate and graduate level students 
to each discipline and, within each discipline, to four student levels on 
the biisis of historical proportions. This was done to reflect the uncer- 
tainty in enrollment projections in the planning process, even though 
actual enrollment data were available for the base year of the runs* 

The complete data set for the validation tests will be referred to 
as the ''reference data set," and the base run using this data set as 
"Run lOA." All sensitivity tests and validation tests involved comparisons 
with Run lOA. The reference data set includes: (1) a 67-'hour week, i.e., 
every hour in which at least one class was scheduled at UCSB in Fall 1969; 

(2) a non-uniform distribution of class hours to the 67 hours, approxi- 
mating the actual distribution at UCSB; (3) the Fall 1969 inventory of 

4 

The methodology consists essentially of multiplying student credit 
hour projections by assignable square feet factors for each discipline and 
space type. The calculations of faculty-generated space are very similar 
to those of the SPACE model, except for the generation of the faculty re- 
quirements themselves. 
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facilities; (4) weekly student hour projections generated by applying a 
single induced course load matrix to the 10-year series of enrollment pro- 
jections; (5) a theoretical class-size-to-room-size assignment matrix;^ 

(6) estimated means and standard deviations of assumed Gaussian dlstribu- 
tions of WSH to class-size ranges, reflecting 8-^ampus histories; and 

(7) non-uniform hour*size factors for lab classes. 

Tables 1 and 2 show the WSH and WRH results for all non-lab classes 
in the validation runs lOA and lOB. The only difference between the two 
sets of data inputs is the starting point for the pseudo-random nuicoer 
generators; all other data are identical. A comparison of the 'actual* 
UCSB data with the outputs of the FAMSIX module shows that, for most of 
the class size ranges, the numbers of WRH correspond closely. in both mag- 
nitude avii variance. For the smallest and largest classes* however, the 
'actual* figures did not include the classes scheduled in locations other 
than classrooms and class labs. The model includes more of the classes 
which are actually taught at UCSB th?.n did the utilization reports which 
were the source of the 'actual* figures. 

It is notable that the average class sizes for the non-lab classes 
vary modestly in the two runs (32,5 and 33.3 in 1969 to 35.7 and 32-9 in 
1978, with 34.0 and 33.3 as the 10-year averages). This variation is 
typical of the historical averages. A similar result holds for lab class 
datta. Also typical of UCSB history (and of the history of the other 

^See pages 11-14 of this report for a discussion of this matrix. 
See pages 6-7 of this report for a discussion of these distributions. 

7 

See pages 8-11 of this report for a discussion of these factors. 
N.B. c =» 1.0 and y « 0 for lab classes; i.e., in the reference data set 
the hour-size factors leave the initial size distributions of lab classes 
unaltered by the scheduling process. 



24 



general campuses) l3 the more marked variation of class sizes within a 
given discipline. Tables 3 and 4 show historical and FAMSIX class 
size distributions in the physical sciences for each of the three course 
levels. 

Comparisons of historical utilization statistics and the utilization 
statistics generated by Run lOA are presented in Tables 5 and 6. It 
should be noted that at least two large rooms which are not classrooms 
are used for large classes at Santa Barbara because the campus has no 
classroom (i.e., no room so classified) with more than 354 stations. Those 
'non-classrooms' are excluded from the 'actual' classroom utilization 
statistics, and thus also from the classroom inventory in the input data 
to the model. SPACE 'constructed' three large classrooms f' . the model- 
generated classes which were larger than the largest classroom in the Fall 
1969 inventory. 

The SPACE utilization statistics reflect, among other factors, the 
class-siLe-to-room-size assignments considered reasonable by the authors, 
and do not represent an attempt to reproduce UCSB's actual utilization 
rates. 

The second type of validation test, comparison of the model's results 
with the traditional U.C. methodology, was performed by excluding the 1969 
inventory from the reference data set. Because of differing assumptions 
regarding enrollment growth a meaningful comparison was possible only for 
the year 1970-71. Only for this year were enrollment inputs nearly iden- 
tical in total for both sets of calculations. Table 7 shows the results 
for the model and the calculations of the U.C. space planners. The siiailar- 
ity of the two sets of results is considered by the authors to be adequate 
evidence of the viability of the model's methodology for augmenting the 
planning process. 
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Sensitivity Tests 

The effects of the randomization procedures within the model on the 
final outputs have already been alluded to in the previous sections. One 
criterion for measuring the acceptability of class size perturbations is 
based on the similarity of the range rf SPACE^induced variations from one 
run to the next with the range of historical variations from one year to 
the next. It might be argued that the actual historical variation in class 
size is an inappropriate base for making comparisons of Monte Carlo per- 
turbations because of enrollment changes within the historical period. Two 

factors mitigate in favor of the comparison however. First, the only prac- 

i 

tical basis of comparison is historical variation. It is obviously impossibl 
to reproduce exactly the conditions which lead to any historical outcome 
preparatory to experimenting with the actual system. Certainly, however, 
for any given student population, the class size distribution within any 
discipline and level is determined stochastically. Secondly, there is no 
clear trend in historical class-size averages induced by enrollment increases 
Presumably any of the historical class sizes are possible for any of the 
historical enrollment totals; and consequently historical variations pro- 
vide a good clue to the possibilities for any single year. 

Tables 1, 2, 3, and 4 display the effects typical of a change 
in the starting point of the randt H'-number generator (Runs lOA and lOB) 
on class size distributions. There is considerably greater variance 
within a discipline and level than there is overall. For non-1 classes 
the overall change in average class size from one run to the other averages 
2.1%. The smallest change for any year is 1%; the largest change is nearly 
6.5%. 



if 
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The ten-year average of class size percentage differences for non- 
lab, graduate, physical science classes is 16%. The smallest change for 
any year is 3.4%; the largest change is 50%. These run-to-run variations 
are generally consistent with the 1967-1969 historical variations for 
this discipline, level, and class type. 

The arguments legitimizing class-size comparisons of historical with 
Monte Carlo variations will not hold for resource needs. Historical data 
on resources expended at Santa Barbara reflect decision makers* prefer- 
ences and political accident as much as actual need. In the test runs 
made so far, SPACE has computed faculty requirements on a workload basis, 
and space requirements of faculty and graduate students on the basis of 
CCHE standards. No effort has been made to simulate the politics which 
generate resources — only the operations which determine resource require- 
ments. Nevertheless, variations of calculated resource needs are an 
important consideration in determining the acceptability of the model for 
planning. Variations much greater than 5% over the 10-year period seem 
unacceptable to the authors. 

Tables 8 and 9 display summary data for Runs lA and IB.^ Faculty 
requirements differ by 1.8%, classroom ASF-years differ by 12.2%, class- 
lab ASF-years by 1.0%, total ASF-years by 2.8%, salaries and support costs 
by 1.8%, total debt-service payments- by 2.4%, and total variable costs by 
2.0%. The 12.2% difference in classroom ASF-years is the only questionable 
variation. Run IB constructed 3 more large classrooms than Run lA; this 
accounts for the difference in floor area requirements. 

g 

Runs lA and IB are identical to lOA and lOB respectively except for 
the exclusion of the 1969 inventory from the input data. These runs show 
what, the model would build to handle UCSB^s enrollments if it were start- 
ing a new campus. A comparison of these runs with the actual inventory 
gives some idea of excess physical capacity i^ the CCHE space standards can 
be assumed to be reasonable* 
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Examining Policy Options 

The usefulness of a simulation model depends in part on the ease with 
which alternative policy what^if questions can be posed and answered. The 
fewer parameters which must be altered, and the easier it is for adminis- 
trators to develop an intuitive feel for what a parameter represents 
operationally, the more useful the model will be, SPACE was initially 
conceived to provide answers for one class of conditional questions, viz: 
What cost consequences would be associated with a change in the hourly 
pattern of class scheduling? The possibility of answering other what-if 
questions arose during the development of the model. 

During the period of development and initial testing it was found that 
it is relatively easy to examine the resource Impacts of adding or deleting 
academic programs, changing faculty workloads, or encountering significant 
shifts in enrollment patterns. In each case it was only necessary to make 
minor changes to one or two sets of data inputs. These data changes 
satisfied the criterion of reflecting operational charges lii an intuitively 
identifiable manner « 

Most of the what-if tests were concerned with changes in the patterns 
of class scheduling. Summary results of three such tests are tabulated 
in Tables 11-13. Table 10 is a summary of Run lOA, in which the 
class scheduling pattern approximated the actual pattern in Fall 1969. 
Changes in the class scheduling pattern were simulated by (1) changing 
the number of hours per week in which classes are scheduled; and (2) chang- 
ing the relative proportions of classes scheduled in those hours • For 
Run lOA a nonuniform distribution of classej to the 67 scheduled hours 
of the week put approximately: 
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1/193 of the classes in each of the 23 least popular hours; 
2/193 of the classes in each of the 6 next least popular hours; 
3/193 of the classes in each of the 8 third least popular hours; 
4/193 of the classes in each of the 16 fourth least popular hours; and 
5/193 of the classes in each of the 14 most popular hours. 

For Run llA, 1/67 of the classes were scheduled in each of the 67 
hours. This *flat' distribution was duplicated for Run 12A but only for 
the 44 most popular hours, i.e., 1/44 of the classes in each of those 
hours. A * lumpy ^ distribution similar to that of Run lOA was used for 
Run 12B, but again, only for the 44 most popular hours, e.g., 2/170 in 
each of the second least popular hours and 5/170 in each of the most 
popular hours. 

The net effects of these appreciable schedule changes were negligible 
in terms of costs. In each case, class sizes declined in comparison to 
Run lOA and operating costs increased correspondingly. Facilities costs 
(debt service) changed only slightly. In general, flat distributions were 
more costly than lumpy distributions and the 67*hour week was more effi- 
cient than the 44-hour week from the standpoint of classroom utilization. 
The analysis of these and other such tests shows that as long as faculty- 
staffing is determined on a contact^hour basis, and as long as students 

9 

show even modest preferences for some hours over others, then even sub- 
stantial alterations in classroom scheduling will do little or nothing to 
change costs. 



Input parameters representing very conservative estimates of the 
class-scheduling effects on class size were included in the reference data 
set and these were left unchanged for the comparison runs. During sensi- 
tivity tests, it was found that still more conservative estimates provided 
little hope for cost-savingB resulting from changing class-scheduling 
patterns. 
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Table 14 summarizes the 10-year projections of six test runs. Note 
that the 2% change in total variable costs induced solely by random varia- 
tions (Runs lA and IB) is greater than the percentage variation induced by 
a scheduling change which cuts by one-third the number of hours in which 
classes are scheduled (Runs lOA and 12B). Care should be taken in the 
comparison of single examples of what-^if tests because of the range of 
results induced by random variations. It is possible that the low end of 
one range would overlap the high end of another. The results displayed in 
Table 14 are, however, typical of the many test runs which were made and 
are representative of the ranges of possible outcomes. 
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FUTURE DEVELOPMENTS 

Planning models, like freeways, are never really completed. There 
are always old sections to be modified and new ones to be added. SPACE 
is no exception. One major weakness of the model is the non^interchange- 
ability of class labs from one level to the nei.t and from discipline to 
discipline. Originally this restriction was intenaod to serve as a surro- 
gate for the actual non-interchangeability of lab facilities between 
fields within disciplines. However > for campuses which regularly schedule 
a certain percentage of their lab classes in class labs of the 'wrong* 
discipline, the current procedure is unnecessarily restrictive. Accordingly 
the Office of Analytic Studies, California State University and Colleges, 
is currently modifying the FAMSIX module to permit user-defined inter^ 
changeability of class labs. 

Another questionable characteristic of the model is its inability to 
construct class facilities optimally. Occasionally the model will encounter 
a need for a new room of, for example, size 9 in one hour and a need for a 
new room of size 10 in the next. Because the model, in its current form 
cannot look ahead at construction time, it would build two rooms in this 
case when one would have served the purpose. By assigning the lowest 
priority numbers, i.e., the highest priority, to the most popular hours 
this problem can generally be: avoided • However, an additional routine to 
'take a second look* at the construction needs should be developed. The 
California State University and Colleges* Office of Analytic Studies is 
working on such a routine, although it appears to be expensive to run. 

The staffing calculations in the VARICOS module are still somewhat 
simplistic. A more realistic model would include consideration of tenure 
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status, death and retirement rates, and rank advancement probabilities. 
These features are still under consideration. 

Finally, at least in principle, the model is adaptable to micro 
analysis of individual disciplines by treating the 33 discipline* cate^ 
gories as individual departments or other course groupings with closely 
related facilities needs. The feasibility of using the model at that 
level of detail will be the next subject of investigation. 
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APPENDIX 



ERIC 



INPUT SPECIFICATIONS: FAMSIX MODULE 



1. IG, IX, SlOl, S210, S400 FORMAT (219, 311) 

Two odd valued integers of 9 or fewer digits, IG is the starting value 
for the randomization of the distribution of weekly student hours to 
class sizes. IX is the starting value for the rotidomization which assigns 
classes of a given size to room sizes. SlOl, S210, S400 are print switches 
for Tables 101, 211 and 212, 401 and 402, respectively. Value = 1 means 
print, value = 0 means don't print. One card. 



2. NCAMP, NRUN, NDATE, NUMYRS, NYEAR FORMAT (4A4, 12, 14) 

Descriptors of the run. NCAMP is a 4-letter abbreviation for the campus 
name. NRUN is a 4-character alphameric code for the name of the run. 
NDATE is an S^haracter alphameric code for the date the run is performed. 
NUMYRS is the number of years for which data is provided in this run. 
NYEAR is a 4'-digit integer representing the base year for the run, eVg. , 
1970. One card. 



2a. DNAME FORMAT (20A4) 

4 cards. The first 3 cards contain 10 entries each, the last card 4 entrie 
Each entry is an 8-character alpha code for the name of a discipline. The 
entries are arranged in natural order, 1-34. The 34th entry must be the 
word "ALL/* The entries should be left-justified in the 8-character fields 



3. ACS FORMAT (15FA.0) 

The average class size in each class size range. One card, the entries of 
which read from left to right in ascending order of sizes. 



ITEMS 4 - 14 PERTAIN TO NON-LAB DATA 



4. NSTATS 

One card with one entry. The number of 
to the number in data set 6). 



FORMAT (13) 
cards in data set 5 (which is equal 



5. WSHR, J, K FORMAT (15F4.3, 6X, 12, II) 

For each combination of discipline, J, and course level, K, one card is 
possible. Each card contains, in addition to the identifiers J and K, up 
to 15 values. These are the means of the distributions of weekly student 
hours to the 15 possible class size ranges in ascending order of size. For 
combinations of J and K for which all values are zero, no card is necessary 
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6. WSHD, J, K FORMAT (15F4.3, 6X, 12, II) 

Same as data set 5 except the entries are standard deviations for the 
distributions rather than the means of those distributions, (The distri- 
butions are assumed to be Gaussian). For each card in data set 5 there 
must be a corresponding card in data set 6, 



7. NOURS, NGRUMP FORMAT (213) 

NOURS is the number of hours in the week during which classes are 
scheduled. The maximum value is 96, NGRUMP is the number of priority 
listings, i.e., the largest no ber in the priority schedule to which an 
hour of the week is assigned. The maximum value is 250. 



8. GRUMP, GlIUMP2 FORMAT (213) 

There are NGRUMP cards, each with 2 entries. GRUMP i.s an integer in the 
interval 1 through 96. GRUMP2 is the priority assigned to the hour (an 
integer in the interval 1 through NGRUMP. The cards are to be offered by 
hour in ascending order. Where there are multiple priority assignments 
to a given hour the cards may be ordered in ascending values of GRUMP2. 
For example, a set may begin with the ordered pairs; (1,6), (1,19), (2,5), 

(2,7), (2,35), (4,16), (5,1), (5,17), All values of GRUMP2 in the 

interval 1<GRUMP2<NGRUMP must be included. No values of GRUMP2 may be 
repeated . 



9. YELP FORMAT (FIO.O) 

The exponent for the function relating class size to the hour of the week 
and the proportion of weekly student hours scheduled in classes during that 
hour. The form of the relationship is: 

X » afp^ . 

Here, x is the average class size in the hour; 

a is the average class size in all hours; 
f is a factor determined for each hour; 

p is the proportion of weekly student hours in scheduled classes 

during the hour; and 
y is an*^ exponent which is not hour-specific. 



10, FACT FORMAT (10F7.6) 

The factor, f, in the function above. There are NOURS such values with 10 
values per card. The first value on the first card is the factor for hour 
#1. The first value on the second card is' the factor fot hour #11, etc. 
Thus, if NOURS « 78 for example, there would be 8 cards. 
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11. DS8, L FORMAT (15F4.3, lOX, 12) 

Proportion of classes of a given class size range which are assigned to rooms 
in each room size range. For each class size range (i.e., 15) there is one 
card. The cards are ordered in ascending values of class size range. L 
is the class size range identifier. The sum of the 15 values on each card 
must be 1.0. 



12. NLREQ, NSTA, NLASF FORMAT (12X, 3112) 

Base year inventory. There are 15 cards, one for each room size range. 
The cards are ordered in ascending order of room size. Each card has 3 
values. The first, NLREQ, is the number of rooms in the size range. The 
second, NSTA, is the total number of stations in the rooms of the size 
range. It is not the number of stations per room. NLASF is the total 
number of assignable square feet in rooms of the size range. 



13. NLFNS FORMAT (1514) 

The number of stations in rooms which will be constructed by the model. 
There is one card with 15 entries. Each entry represents the stations per 
room for a given room size. The entries are ordered from left to right in 
ascending order of size. 



14. FFPSNL FORMAT (1514) 

The number of assignable square feet per station for rooms to be constructed 
by the model. There is one card with 15 entries. Each entry represents 
the asf per station for a given room size. The entries are ordered from 
left to right in ascending order of size. 



15. 22. 

These data sets are identical respectively to data sets 4 11 except that 
they represent LAB classes and rooms. Note that NSTATS, NOURS, and NGRUMP 
may have different values for lab and non-lab data. 



23. NCARDS FORMAT (13) 

The number of cards in data set 24. 



24. L, LREQL, NSTAL, LASF, J, K FORMAT (12, lOX, 3112, 20X, 212) 

For each class-lab size, L, discipline, J, and level, K, for which there is 
non-zero data, there is one card. LREQL is the number of class-labs in the 
base year inventory for the L-J-K combination. NSTAL is the total number 
of stations in those class-labs* LASF is the total nuniber of assignable 
square feet in those class-lahs. The order of the cards is not important. 
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25. NCARDS FORMAT (13) 

The number of cards in data set 26 (which is equal to the number of set 27). 

26. LFNS, J, K FORMAT (15lA , 8X, 212) 

For each combination of discipline, J, and level, K, with non-7.ero data 
there is one card. Each card contains the number of stations in class-labs 
to be constructed by the model in each size range. There are 15 values per 
card, each representing the stations per lab for the 15 size ranges ordered 
from left to right in ascending order of size. The crder of the cards is 
not important. 

27. FFPSL, J, K FORMAT (15lA, 8X, 212) 

For each card in data set 26, there is a corresponding card representing the 
assignable square feet per station in labs to be constructed by the model. 
The sets are identical in form and specification except asf ^per-station is 
►to replace stations=per-lab.' - 

ITEMS 28 AND 29 ARE TO BE REPEATED FOR EACH YEAR OF THE RUN (i.e., 28, 29, 

28. 29, 28, ) 

28. WSH, J FORMAT (3F10.2, 39X, 12) 

Non-lab weekly student hours. There are 33 cards, one for each discipline. 
If there are no weekly student hours, for a given discipline a card which 
is blank except for the discipline number in columns 70-71 must be included. 
The cards are ordered in ascending order of discipline. Each card contains 
3 values, one for each course level. The first value is for level 1, the 
second for level 2, the third for level 3. 

29. WSH, J FORMAT (3F10.2, 39X, 12) 

Lab weekly student hours. This set is identical in form and specification 
to set 28 except Non-lab is replaced by Lab. 



NOTE: The maximum number of disciplines is 33. 
The maximum number of levels is 3. 
The maximum number of sizes is 15. 
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i.N^UT SPECIFICATIONS: VARICOS MODULE 



It is assumed that the following shall apply to card-input. If VARICOS is 
run in tandem with FAMSIX, the user may ^'equest that FAI4SIX will write data 
sets 1, 13, and 14 onto tape or temporary disk storage. If this is done, 
it is assumed Chat the storage file will "look like cards." 



1, NCAMP, NRUN, NDATE, NUMYRS, NYEAR FORMAT (4A4, 12, I4) 

1 card, NCAMP is a 4-character alpha abbreviation for the campus n^me, 
NRUN is a 4-character alpha abbreviation for the name of the run. NDATE 
is an 8-character alpha abbreviation for the date of the run, NUMYRS is 
the number of years for which data is provided. (Maximum value is 10). 
NYEAR is the base year for the run (e,g,, 1969), 



2, DNAME FORMAT (20A4) 

4 cards. The first 3 cards contain 10 entries each, the last card 3 entries. 
Each entry is an 8-character alpha code for the name of a discipline. The 
entries are arranged in natural order, 1-33. 



3. Al, A2, Bis FORMAT (6F10.0) 

132 cards arranged in 33 sets of 4 cards each. The first card contains Al, 
•the second/A2, the third Bl, the fourth B2 for each set. The sets are 
arranged bW discipline in natural order. Each card contains 6 values: 
first Non-labvdata for course levels 1-3 followed by Lab data for course 
levels 1-3. Al is the proportion of weekly room hours in the appropriate 
level, type, and discipline which are taught by TA's, A2 is the same 
data for the non-TA part of the academic staff. (If there were no team 
teaching -il + A2 would equal 1.0) < Bl is weekly contact hours per FTE 
for TA's, B2 is the same data for the non-TA part of the academic staft. 



4, FACSAL FORMAT (6F10,0) 

1 card containing 4 values. Average yearly faculty salaries for the four 
faculty types. In ord*?r, these are: TA's, irregular ranks, regular non- 
tenured , tenured . 



5. C2, D2 FORMAT (6F10.0) 

33 cards arranged in natural order by discipline. Each card has 6 values: 
2 for C2, and 4 for D2. The first C2 value is non-TA instructional staff 
per TA in non-^Lab courses, i,e,, supervision by non-TA faculty of TA's 
generated by non-Lab classes* The second C2 value is the Lab--class counter- 
part of the first, (Typical values are 0,1), The D2 values are non-TA 
faculty per student for each student level: LD, UD, Grad 1, Grad 2 in order. 
These values represent independent study workloads when A2, B2, and C2 are 
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5. (continued) 

non-zero. They may be used to generate faculty on a strict student/ faculty 
ratio by setting A2, B2, and C2 to zero. However, it is necessary to 
generate TA's via workload computation and the Al, Bl parameters, 

6. SUPRAT, RKDIST FORMAT (6F10,0) 

33 cards arranged in natural order by discipline, 4 values per card. The 
first value is discipline support $ per non^TA faculty. The next 3 values 
are the distribution of non-TA faculty by rank, i,e,, proportion of non-TA 
faculty which are, in order. Irregular, Non-tenured, Tenured, 

7. Rl, R2, 01, 02, SI, S2 FORMAT (10F8,0) 

33 cards arranged in natural order by discipline, 10 values per card. The 
first 2 values are research ASF/non-TA faculty and research ASF/TA, Then 
research ASF/Grad 1, and research ASF/Grad 2, 01 and 02 follow the same 
pattern replacing research ASF \?±tb office ASF, The ninth value, SI, is 
support ASF/research and office total ASF, The tenth value is* support ASF/ 
class-lab total. 



8, CAPFAC FORMAT (10F8,0) 

33 cards arranged in natural order by discipline, 5 values per card. 
These are capital cost .factors for each tirpe of space, i,e,, $/ASF for 
classrooms, class labs, ;esearch, office, and supporting ASF, Each of 
these numbers will be multiplied by an index number (cf , Item 11 below) so 
it is possible to treat these numbers as multipliers of an overall cost/ 
ASF index. 



9, K, INASF (3-5, K) , NDCHK FORMAT (12, 3110, 12) 

Inventory of non-class ASF, For each discipline, K, with non-zero inventory 
of non-Class ASF there is one card, !^»ach card contains in order: discipline 
number, inventory research ASF, inventory office ASF, inventory support ASF. 
The last physical card must additionally contain a 1 in cc, 34 to signal 
the end of the set. 



10, SINFL, BR, BL FORMAT (2F5,0, 12) 

One card containing 3 values* SINFL is an inflation factor - annual rate of 
inflation for physical facilities less annual rate for operating costs plus 
one - e.g,, 1,01, BR is the rate of interest on bonds, e,g,, 0,03. is 
the bond life> e*g., 25* 



11, CSTDST, D4MN0, DPSQFT FORMAT (10F8,0) 



ERLC 



One card containing 6 values. The first 4 values are the aistribution of costs 
over the 4-year period of capital construction for any project, i,e,, the 
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11. (continued) 

proportion of the total project costs incurred in the first through fourth 
years of construction. The fifth value is the overall $/ASF figure for 
maintenance and operation of physical plant. The last value is the overall 
$/ASF figure for construction (or the comparable index figure - see Item 8 
above) • 



ITEMS 12-14 ARE TO BE REPEATED IN ORDER FOR EACH YEAR OF THE TEST RUN 



12, ENR, NDCHK, K FORMAT (AF5.0, 12, 46X, 12) 

Enrollments by discipline and student level. For each discipline with 
non-zero enrollments, one card containing 5 values. The first four values 
are the enrollments for LD, UD, Grad 1, Grad 2, in order. The fifth value 
is the discipline number (cc. 69, 70). The last physical card must contain 
a 1 in cc. 22 to signal the end of the data group. 



13. I, K FORMAT (212) 

WRH FORMAT (3lA) 

Weekly room hours by course level, discipline, and class type. For each 
combination of discipline and class type with non-zero weekly room hour 
data there are two cards. The first card contains two values, the class 
type, I, and the discipline number, K. The second card contains three 
values: the number of weekly room hours in each course level for that 
discipline and class type. The WRH are ordered LD, UD, Grad. 

At the end of this group of cards, there must be one blank card to signal 
the end of the data group. 



14. K, INASF, NWASF iORMAT (12, 4110) 

33 cards containing 5 values per card. The cards are ordered serially by 
discipline. The first value is the discipline number. The next two values 
are the non^Lab and Lab inventories. The final two values are the non-Lab 
and Lab newly constructed ASF. The inventories include the new ASF. 
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FAMSIX - PROGRAM LISTING ~ 



c 

C J HIS IS FAMSIX OESIGhED AND WRITTEN BY DQM^M SMTHt W. GARY WAGNliK, 

C hNO JOHN LAFLtR OF THh CFFlCf: dF ANALYflJAC bTJUltSt 

L VICE PKESIUENT - J^LANNINlit UNIVERSITY UF LALlFJKi^lA. 
C 

C THE CODING OF SUBRUUTINE NLSPiL HAS b£EN ACAPftD F8Jf1 CCH£~FAM» 

C A PRUDUCT UF MATi iEN AT ICA INC. 

C 

C COMMENTS REFER Jhll FLOW CHART OF FAMSIA dV BOX Lt:TTER* 
C 

C THIS PROGRAM HAS UEfcN TIMED BY A 60GLE A^0 ijAtibAGE PkUOUCT AMO 

C FUt-TUNir4G PROGRAM CHANGES HADE FOR INCREASED EFFICIENCY. 
C 

C WEUKLY STUDENT HOURSt CLASS SUES AND ftCGM HOURS AKb GENERALLY 

C CARRIED AS INTEGERS ^ AO ROUNUED WHEN NECESSARY. 

C 

C INTEG£R*2 IS USED THROUGHOUT THE PROGRAM TO MINIMIZE STORAGE* 
C 

C HEkE FOLLOWS A DICTIONARY OF HOST VARIABLES USED IN THIS PROGRAM 

C tXCGPT THOSE OF THE FORM L,LLfJJJ ETC WHICH AR^ LOCAL COUNTERS 

C uR DO LOUP ARGUMENTS* 
C 

C ACS IS Trit AVERAGE CLASS SIZE FOR 15 CLASS SIZE INTERVALS 

C ADUM IS THE TEMPuRARY HOLDING VARIAGLii HOLUiNii GAUSSIAN RANDOM NUMBER 

C AFT THE WORD •AFTER* FOR TITLES 

C AWSH IS THE INITIAL COUNT UF WEEKLY STUDENT HJURS 

C BEF THE WORD •BEFORE' FUR TITLES 

C DATA TEiMPUrtAHY VALUES ACCUMULATED FOR P^IINT LlNl:S 

C CNAME OiSClPLlNfc NAME FOR REPORT PRINT Ii>iG 

C CSii (FAM DATA SET Q) SHOWS PRubAolLlTY JF CLASS UF A GIVEN SIZE 

C ciiiING IN ROOM OF A GIVEN SIZE 

C FFPSL SP \CE standard: NUMBER ASF PER STATION rOK LA3S TO BE BUILT 

C FFPSNL SPACE STANDARD: NUl^rtER ASF PER SfATiJN i=OR CL/SS RU,OMS TO 

C aE BUILT / 

C FPS TEMPORARY VARIABLE: FEET PER STATlQ;^ USED IN TAB LES /30i , 302 

C FPWSH ASF PER WEEKLY STUDENT HOUR / 

C GASP IN HDU:^ LUCP IS TEMPORARY VARIABLE HOLDING DISTRIBUTION OF 

C CLASS Sl/.ES BY ROOM SIZES 

C IBLD TEMPURARY NUMBER OF ROOMS TO BE liUlLT IN HJUR LOOP 

C IG IS THE VALVE USED IN GENERATING RANuO.-l NUMtitaS FOR THE 

C VARIATION JN WEEKLY STUL>ENT HOUKS. 

C INASF CURRENT NUMBER OF ASSIGNABLE SQUARE FEET FOR USE BY VARICOS 

C IX IS THE VALUE USED IN GEULRaTING RANDOM NUMBERS FOR THE 

C VARIATION IN ROOM SIZE ASSIJNLD TO A JIVEN CLASS 

C FACT IS HOURLY FACTOR USED IN ADJUSTING CLASS SIZE BY HOUR 

C GRUrIP IS THE ■HOUf< TO WHIUI EACH PRlUkXTY IS ASSIGNED 

C GRU:*^P2 IS THE PiUORITY ORDER OF THIS i>RIORlTY 

C JUM IS A TEMPORARY VARlAtiLt i<t PRESENT I NG CELLS UF JWRH IN TURN 

C J WALK IS AN ORDtKLiD LIST OF CLASS HOURS 



JWRH IS USED TWICE. THE FIRST TIME IT IS THE .^£cKLY ROOM HOURS 
UER WED . FROM THE ORIGINAL WEEKLY STUJENT HJO.<S, THE SECOND 
TIME IT IS THE WEEKLY ROO.M HOURS ASSIGNED tN THE HOUR SIZE 
FUNCTION BY DISCiPLlNEf LE/ELt AND SIZE FJR fHAT HOUR. 
KHCSZ SUMS OF CLASSES BY CLASS SIZE FOR TALLtS ^0i,40Z 
KHRSZ SUMS OF CLASSES BY RUUri SIZE FOR TABLES 40i,H02 
KWRtI HOURLY CLASS HOURS BY .01 S C I ?t INE , LEVEL AND REVERSE SIZE 

ALL CAREFULLY LUMPED AS SINGLE. SUBSCRIPT 
KKWUH HOURLiT CLASS HOURS BY SIZE ONLY, JSED ^^^ HOUR SIZE FUNCTION 
KWRHD SAME AS KWRH OUT PKOPERlY SUtiSCKIPTEU 



C 

C 
C 
C 

C 

c 
c 
c 

c 
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C LASF NUMbfcrt OF ASSIGNABLl: SQAkG ft£T IN LABS 

C L8LUL LAO RUGMS TO aE BuILT 

C LFN3 IS THE NUMbER OF STATIONS SUt FGK CLASS LABS TO BE 

C 3U1LT 

C LP IS USiiO FIRST AS THt SIMPLIFIED SUtiSCKl P T , JF FIRST JWRH, LATcR 

C IT IS POINTilK TO CURRhNT PKiU^lTY T 0 Bt USUJ IN HOUR LOOP, 

C LPP IS UOuO AS POINTER TO EVEKY CLASS HOOK WHctN MAKING UP JWALK. (LP 

C IS ONLY PQIMfhR TO EVERY CELL OF CLASS HOURS,) 

C LREQ L^SF SUhHLD OVER DISCIPLINES FOR PRINTING 

C LREUL iS THE INVENTORY OF CLASS LAOS AT ANY TIME 

C MREQ IS LREJ)L SOMMEO OVER DISCIPLINE FOK RtrMRT r^URPOSES 

C MWRH CLASSES OVER ALL MOORS dY UiSClPLlNEt LEVEL, CLASS SIZE 

C .ANQ ROOi^ SIZE FOR LABS, 

C HhRHt) TGMPCRAKY COUNT OF ROOKS 6Y SIZE, Uf^ClPLlNE AND LEVEL 

C WITHIN HOUR LCOP, UPOATJfU IN NLSPIL. 

C M Is THE NUMBER 1 REPRESf:NT INO THE CARO REAi;GR IN INPUT 

C N5 IS Trit NUMBER I Kef' RES ENT INO THE OUTPUT FILE FUR VARiCOS 

C N6 IS THE NUMBER 6 REPRESENTING THE PRI.^iTER lA OUTPUT 

C NCAHP IS NUMdER UF THE CAMPUS (FOR PRINTING} 

C NUATA TE*S. JRARY VALUES ACCUMULATEU FOK PRINT LJNES 

C NCATE IS THE OATc(a CHAR FOR PKINTING) 

C .NERR IS TiiE NUMQER OF ERRORS FUUNO DURING READ IN OF OAT A 

C NEWRH TOTAL \^EEKLY ROOM HOURS OVER ALL dOJRS. 

C NEViTOT IS FIRST JV.'RH SUMMEO OVtK DISCIPLINES FOR REPORTS 

C NGRUMP IS THE NUMBER OF PRIORllItS IN OCCURENCE FOR hOURS 

C NIT IS THE ITERATION NUMBER (FUR PRINTING) 

C ^LASF TEMPORARY VARIABLE: USEO FOK CCi-IPUTlNG UTILIZATION RATES IN 

C TAaLES 301,302 

C iNLBLD NON-LAB ROOMS TO BE BUILT 

C NLFNS IS THE NUMBER OF STATIONS BY SUE FOR CLASSROOMS TO 8E 

C BUILT 

C NLRLU IS THE INVENTORY OF CLASSROOMS AT ANY Tl«t 

C NNHUlO NUMBER OF STATIONS: URiGINAL + BUILT BY SIZE AND LEVEL 

C NOl TABLE NUMBERS FOR TITLES 

C NONK.H SUM OF NEWRH OVER DISCIPLINE FOR TOTAL LINES 

C NOUKS IS THE NUMtJER OF HiJUHS WHEN CLASSES HAY OCCUR 

C NPHU "NON* OR BLANK FOR TITLES 

C NRMS THE SUM OF THE WEEKLY ROOM HOURS li^ THAT HJURiONE NUMBER) 

C NRCCM •ROOM* OR *LAB« FOR TITLES 

C NRUN IS THE RUN NUMBER ( H OR PRINTING) 

C NSTA IS THE NUMi3bR OF STATIONS BY SIZE FOR CLASSKOOMS 

C NSTAL IS THE NUKbER OF STATIONS BY SIZE FOK CLASS LABS 

C NWASF ASSIGNABLE SOUARt FuET JF NLW RLCMS L LA-iS FOR VARICOS US£ 

C NYEAR IS THE NUH'icR OF THE YEAR (FOR PRINTI<>iG) 

C SIOl INPUT Sv^ITCH TO PRINT TABLE lOl 

C S210 INPUT Sv^iTCH TO PRINT TAbLbS ^.11,212 

C S^OO IfyPUT SVitTCH TO PRINT TA.H.LS ^01,402 

C WRHS TCNPOKARY COUNT Of KUOMS BY SIZE WITHIN HOUR LOOP, DPOAXEO 

C IN r.LSPlL. 

C V^RHSAV IS HEtKLY ROOM HOURS ASSIGNED IN THt HOUR SIZE FUNCTION BY 

C UlSCIPLINEt LEVEL AND SIZE FOK ALL HOURS 

C hRhSUM TOTAL Vt'EEKLY ROOM HOUR!^ OVER ALL HOURS 

C WRHZUH CLASSES OVER ALL HOURS taY CLASS SIZE AMO ROOM SIZE FOR 

C NUN-LAliS, 

C hSM IS THE BASIC WEEKLY SfUUENT HOURS BY, 0 1 SC IPL I NE ,LEV£L 

C ANO LAB-NUN LAB AS READ IN <NOT OESfROYiiU) 

C HSML IS THE PANOuMIZED WSM INPUT OVER CLASS SIZE, 

C . ^ShO IS THE DEVIATION USED IN GENERATING RANOJM CHANGES IN 

C WEEKLY STUDtNT HOOKS 
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C ViSHR IS THE MtAN Of THE OEVIATJUN USED IN Ge«4SKATING RANDOM 

C CHANGtS IN WLtKLY STUOfcNT HOURS 

C XDUM dtCOMES THE CGKRECTiUN KATIU TO STAaELlZc THG WEEKLY 

C STUUEiU HOUKS liY OISCXHLlNfc ANU LEVcL IN SP I ft OF RANOUM- 

C UING CLASS Slit: 

C YELP IS THE POWER FACTOR UStO IN ADJUSTING CuASS SIZE 5Y HOUR 
C 

UIHc^SlON ACSU5)f DSO 1 1 5 , 15 , 2 ) , NOATEl^)f HSH(33f3,2li 

1 AIJUM(15}, F ACT ( 96,2) , NLASFl 15) , WS.iU( 3 3,Jfl5,2) i 

2 AFT(2) , INASFl2,33i, NPR012), «ISHK( 33,3tl5, 2) t 

3 AWSH(3J,i,21 fLASFI 15 1 33, 3 ) i NRCCH ( 2 ) , WSH 1 ( 33 ,3 1 1 6 1 » 
^ BEF(2), LREQ13,15), NWASfM 2 , 33 ) , YtLH ( 2 ) , 

i ^ 0ATA{3}f N0ATA(3), NOKoJi 

6 DMAM£(2f 34) 

C 

iNTEGER*2 FFPSLC 15 #33, 3) t KWKHO ( 33 t 3 t 1 5 ) , NNHOLO (3 r 15) , 

1 FFPSNLCL5) . LBLOLl iS,33#31 , NONK.H 3,15) , 

2 GASP115, 15) , LFNS< 15,33. 3), NOUKSC 2) , 

3 GRUMP(250,2) , LREQL 1 15 , 33,3 ) , NSTA(15), 

4 GKUMP2 (250,2) f MKhyi3,15), NSTALdS ,33 ,3 ) , 

5 IdLD{ 15) , MWKHdS, 15f 33,3J , k.^HSUSIt 

b hWRHO( 15,33,3) , rtilHS AV ( 3 , L5 , 33 J , 

7 JWALKl 12500) , NE WRH ( 33 , 3 , 15 ) , WKHS0iiU5) , 

d JWHH(3, 15,33) , NEWTOTU, 15) , WRHZJ.-I 1 15 , 15 ) , 

9 KHCSZU5) , NGRUMPiZ) , 

A KHRSZ(i5), NLBL0115), SLOlt 

d KKWkHll5), NLFNS(i5), S210, 

C KWRH(I48 5I , NLREQ( 151 , 5400 

C 

C EQUIVALENCE STATEMENTS ARE FOR oTORAGE EFFICIENCY EXCEPT FOR VARIABLES 

C KWRH ANC KWftHO WHERE A FUNCTIONAL CROSS rXlSTS. 

C 



bUUI VALENCE i AOUM( 1 ) ■ NEW RH { 1 ) ,WRHSUH(1J ) 
EQUIVALENCE <KWRH( 1) iKWKHOt 1) ,MWKH0(1) ,WRHS( 1) } 
EGU1VALENCE(MWRH(U,WRHZUH<1) ) 
EQUIVALENCE (MW^H(6001) ,WSH1(1) ) 
cQUlVALENCE C LtiLDL U ) , NLBLO ( 1) ) 

C 

C JATA STATEMENTS ARE FOR REPORT HEADINGS 

a ATA N0I4 1),N01(2),N0i(3},N01(4)tN0i(5} ,N0 1 ( 6) , NKOOM ( 1 ) , NRQOM( 2 ) / 
I'Ol^ ,'02 S 'II* 12» ,'21* I '22* ,'.^COM* , ' LAB'/ 
UATA NPRGl 1) ,NPR0(21/«N0N» ,* '/ 

JATA UEF(i},bEF(2),AFTU)f AFTl2)/"bEFS»0KE' .'AFf ,'ERV 

Nl^i 

N6=6 

W5=^2 

NZERC=0 

NERR=»0 

C 

C FOR SIGNIFICANCE OF INPUT VARIAtiLES SEE INPUT DtJCUMfNTAT ION (BOX A) 

REAU (Nl.^OO) IG,IX,S10l,S2lO,S400 

wRireCNcooi) iG, IX 

READ tNl,902) NCAHP,NKUN, NOAT E ,NUMYRS, NYEAR 
WRIT£(N5,902 )NCAMP,NRUN,N0ATE,NUHYRS,NY6AR 
KEAuiNl,929) DNAME 
Kt/lUlM,J03) ACS 

C 

C - FCLtAR, KCLEAR( INTEGER«2 ) I ANO KKULEAR ARL SIMPLE CLEARING ROUTINES 
CALL fCLtAR( WSHK,2970) 
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CAU FCLEAR( rt5HDf2970) 

CALL KCLeAW(Lj^EQLil485) 
CALu KCLtARlNSTALf 1485) 
uAlu KKCLER(tASF,l463) 
CALL KCLEARCFFPSLf 1485) 
UG 10 n = l,2 
KEAU (Ml V04)NSTATS 

ua i I M=iiNsrArs 

ii kEAJ(N1,^05) JiKi {WSHK<J .K.Lf II I iL'^iff 15) 

JO 12 f<=lfNSTATS 
1^ R;EAU(Nlt^05) JrKi ( WSHC)(JtKffLtU) ft^ltlS) 

KeAJlM,9J4) iNOURSdU fNGRUMP(U) .j 

UGK^NOHUMP a I I 
OQ 13 J^lfNGK 

KtAU (Nlf'^O^) GRUMPC J,lli Gi^UMP2« Jfin 
KEA0(Nlf906) YELPdll 
NGU=NOURS( 11 ) 

KEAj<Nlf907) t FACT IN, I 1) •N^l.NOU) 
UO 1^ LL=l fl5 

Kt^AuCN 1,^08) Lf (0S6(L,M,in fM=l,15) 

C 

C VAIUCUS OATA CHECKS ARE PERFORMED BUT THE P/^OGKAM ONLY STOPS AT 

C titGINNING OP ACTUAL COMPUTATION IP ERROrtS HAVE d^EH FOUND. 

IFtL.NE.LL) CALL ERROR INfcKR, 6 ) 

C 

C bSQ IS MADE INTO A CUMULATIVE LIST 

ua i4 M = 2f lt> 

14 USOIL ,M,II )-OSQtL,M, 11 )+0S8 ( L iM-l • I 1) 
iF{NbRR.N£:.Q) STOP 

IFU U£g.2 ) GO TO 16 
UU 15 L=U15 

15 KEAU<Ni,'?09) NLKEQI L) , NST A I L ) i NLASF ( L J 
KEAu(M,910) NLFNS 

K£AU(Ni,910) FFPSNL 

00 TO 10 
lo KeAU(M,SQ4) NCAROS 

oa 17 M=l, NCAROS 
17 KEAJ(Nl»9in LiJ,K,LREuiL(L,J,K} ,NSTAL(L,JiK),LASF(LiJfK] 

READ (Nl,904J NCAROS 

UO xe H=lf NCAROS 
Id KGAulNl,912) J»K,(LFNS{L,JiK),L=lf 15) 

JQli M=»lfNCAROS 
i.9 K£AU(NI,912) J t K , ( FFPSL ( L , J ,K) i L=l t 15 ) 
10 uCM INUE 

C 

C HERU STARTS MAJOR LOOP OF PROGRAM DOING ALL WORK PQiK EACH YEAR(BGX B) 

uO LOG NIT^UNUHYRS 
IF(NIT.GT» IJ NYCAR^NYEAR ♦ 1 
00 101 II==lf2 
UG IGl JJ-it33 

READ (Nl« 9213) J , ( WSH( J»;<f I 1) >K^ly3) 
IFtJJ.NE.J) CALL ERKORlNERRfS) 
lUi tONi INUE 

U-{.WERR.N6.0) STOP 
GALL rCLEAR( AViSH>19<J) 
CALL KKCLtRi INASF,66) 
CALL KKCL£R<NWASF»66) 
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C TABLE 101 

1F(:>1C1.£Q.0) GO TO 102 

CALL PtlZPKNOlll) fNCAMP,NY6AR,NITf NftUNfNOATCWSHt IfONAME ) 
i^J^ oGNTINUE 

C 

C HERc STA»<TS LCGP FUR CLASS TYPE: Z«LAbf l« NUN-LAd. (BOX 0) 

OQ ^0Q U^li2 
CALL KCLEARi WRHSAVf 1485) 
WoR-'NGKUMP (I [ ) 

U3 ZOl J'l f33 ^ 

00 ^01 K=*lf3 

aDUM=0. 

00 ^02. L^l fl5 

aOUM(L)«0. 

C 

C GAU:>5 IS RANOOM NUMBER GENERATOR FUR NORMAL OBVIATES 

lF(hSHiJ(J,K»LtiI).l^T. .OOOl.OR. WSHR( J t K , L f I I ) • G T • .0001) 
LCALl. CAUSS(IG,WSH0( JfKfLf I I ) t WSHR IJ f K f L , 1 1} , ADUMt L M 
iF(AOUM(L) .LT» 0.) AOUM(L)=0. 
ZOZ aDL'J^ = AOUM(L)+X0UM 

iF(XCUM#GT. .001) X0UM*1*0/X0UM 
UU ^01 L«lfl5 

C 

C UISTKIBUTE WSH TO CLASS SUES BY DISCIPLINE AND LliVEL (dOX E) 

^Oi WSHK J ,Kf L )«A0UM(L)«X0UM«WSH( Jf Kf II ) 
JO 210 J=lf33 
UO 210 K^l ,3 
00 210 L^ltl5 

C 

C HGRHALIZE ViSH AND GENERATE INITIAL WEEKLY CLASS HlURS BY OIVIOING 

C 6Y AVERAGE CLASS SfZL AND ROUNOING TO INTEGERS (BOX F) 

21J JWRH(K,Lf J )^WSHl(JtKtL)/ACS(LI+.5 

iF(S210,Eg.O) GO TO 213 

CALL KCLEAR( NE^T0T,45) 

WHirE(N6,9I3) NOK II+2i,NCAU',NYEAR»NPR0(ll) ,NI T ,j^RUN t NOATEt BEF 
JO 211 J^li33 
JO 212 K^lfi 
00 212 L=lil5 
212 NEWTOT(KtL)=NEttTOT(K,L)+JWRH(KfL, J) 

C 

C TABLES "211,212 

III CALL PEtP2(J»JWRH(l#lf J) fACSfllfOfONAME) 
CALL PEEP2(34,NEWT0TtACSf I I tOi ONAME) 

C 

C THE 00 220 LOOPS AkE IN PREPAKATICN FOR PRIORITY SCHEDULING 

C v^ITHIN TMG HOURLY LOOP. THE VARIABLE KHRH IS THE VARIABLE 

C KWKHO CUNSIOGRtO IN THE FORTRAN STORAGE JKOER AS A ONE-DIMENSIONED 

C VARiADLt:. JWALK IS AN ORDERED LIST OF POINTcRS Tj3 THIS VARIABLE 

C :>CRrED BY REVERSE SUE OF CLASS, LEVEL AND aiSCJPLlNL. WHEN fHE 

C PRIuRITY ASSIGNMENT IS WANTED IT WILL BE EASY TO •WALK* DOWN THE 

C J WALK LIST. 

21J LP=J 
LPP==0 
SUM«0. 

UO ^20 LL^1,15 

00 220 K^lfJ 

UO 4:20 J=lf33 

JUW-=JV^RH(Kfl6-LL,J) 

lP=lP+1 
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lF(JUH.£g.O) GO TO 220 
:>UW=SUH*JUM*ACS( 16-LL) 

^.iJ CUNf INUE 

lF(i I .EQ.^ ) GO TO 230 

lALl KCLtARlhJLdLO, 15) 

CALL KCLEARl WKHZUH,225) 

i-ALl KCLEAKUKHSUMf 15) 

oC TO 231 
2jO oALl kcl^.ar< LBLOL, i^»a5) 

CALL KCLeAKlMWftHf 22275) 

i-ALL KCLEAR(New.^H,l485) 
2ii iF(:^4CQ.&Q.l ) Wi^I Tt: (N6 «9 14IN01 ( I 1 ) •WROUMii t) 

NCU=NGURS( II ) 

CALL KCLEARt IBLU.15) 

CALL KCLcARC JWRHfl485l 

CALL KCL6AR(KWRH, l^a:>} 

CALL KCL6AR(GASP,225) 

C 

C MERt STARTS THE HOURLY LOOP (60X G) 

UO jOO N=lfNQU 

C 

C CHECK TO SEE IF ALL REMAINING HuUKS ARE EMPTY 

1F(lP.GT*NGHUHP(IJ )) GO TO 300" 
NSW=0 

c 

C CHECK TcJ ^:EE if THIS HOUR HAS A PRIORITY ANJ SET >HiTCH 

3v)i iF(^RUMP(LPf 11 ).NE,N)GO TO 303 
NSW=1 

rt = GKUMP2lLPf I I ) 
lP-uP*1 

UC j02 ^^ = ^|LPF^NGR 
KNN^i JWALK(NN) 

C 

C ASSIGN CLASS HOURS BY PR10RITY<S£E LONG COMMENT AJilVt) (BOX H) 

30^ iS.HRH(NNN) = KWRH(NNN) + l 

GO fC 301 
JU3 CONf IWUE 

iF(i>4SH.E0.0} GO TO 300 
iF{£^a0.tQ.O) GO TO 304 
NRMS^O 

CALL KCLEAR(KHCS2f 15) 
CALL KCLLAR<KHRSi,13) 
30-^r »<HSUM = 0 

UC ilO L^i»l5 
NKlir<H(L)=0 
□ 0 j 10 J^lf33 
00 310 K=l,3 

IF IKl^RHUi JfK,16-L) .EQ. 0 1 GO TO 310 

C 

C FINO NUM9i=H UP CLAiiS HOURS BV^ SIZE AND HOUR IN PREPARATION FOR THE 

C HOUR SIZE F-UNCriUN ROUTINE. 

KKMKH(L) ^KKl-JRMt L)+KWKHD( Ji K,l6-t) 

C 

C rtHSJN IS SUM UF KSH FOR HOUR ANO IS NEEUCD TOR FX. 

WHSUM«=HHSUM+KWRHO( JfKf 16--L)*ACS(L) 
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ilU CGNflNUE 

C 

C JHlb G£N£K/iTHS THE HOURLY SUE VAKlAlION 

hX^FACKNf ii J ♦(WH5UM/SUMi'J'<'Y£LP(ll i 

C ]HE hCOH SlZt f^JNCTlON APPLIES FX Ai\U CHAwGiiS kUU.i HUURS CBOX 1) 

CALL HKSZFN(KKWRHf ACSi JWrtHf KrtKMUf IXp lYf FX) 

KSVsCH=a 

c 

C NOTc THAT THIS LONG 400 LOOP IS t3EING DONE BY YEAKf CLASS TYPEf HOURt 

C UlSClPLlNEi AND LEVEL ANO THAT IT HAS UGJbLE ANU TRIPLE LOOPS 

C rilTHIN IT. EVEKY OPPQKrUNITY IS TAKEN TO BYPASSf BUT 

C COOINc; AOOEO HERE CCULO BE EX^»E^^S^VE. 

00 'VOO K^Il ,3 
uG 410 L^ifLS 

if- lJXRH(KtL»J) .EQ. 0) GO TO 410 
NN^^JXRHiKtL, J) 

HRHSAVI K ,L i J) -WRHSAVC K ,L f J ) +NN 
JWRrt (KfLi J)=0 

UO 411 l»i »NN 

C 

C KANJU IS A RANOGM NUMOER GENERATOR FOR A UfUFJKHLY JISTRIOUTED 

C KANOUM VARIABLE ON THE CLOSED 0»i INTERVAL 

v-ALL RANOO (IX, lYfRATIU) 

IX=XY 

□ 0 412 M^L ,15 

C 

C CLASSES ARE ASSIGNED TO ROOMS ACCORDING TO fHc INPUT PR0BA6ILITY 

C OlSTRIBUr lUN (OSS). AT THIS POiNT WE SWITCH PROM WLEKLY CLASS HGUkS 

C rO WEEKLY RCOM HOURS. 

il^(KATlG*GT.OSa(L,M,I U ) GO TO 41^ 

GASP(M,L)=GASP(N,L>+1 

IFCl I.EU.2 ) GO TO 413 

rtRH:>(Mi=MKHS<M)>l 

GO rO 411 
^.^3 rtWPHD(M,J, K)=MWRHD(Mf JfK) + l 

oO 10 ^11 
41^ CCNT INUE 

WRITE(N6,915) 

;>TOP 
411 CONTINUE 
41U CENTINUE 

if ii I .EQ.l ) GO TO 400 

irCKSWCH.EQ.O) GO TO 400 

KSWCH=Q 

c 

C I^LSPIL ASSIGNS CLASSROOMS. ROOMS ARE NOT •liUiLT' IF A LARGER ROOM 

C IS FREEt BUT THE CLASS IS SPILLED UP. (BOX K) 

CALL NLSPIL( lOLDt 1!>,HWRHJ( If JfK) , LREuL ( 1 « J tK ) , GASP ) 

C 

C LGLlECT DATA FOR TABLE 401 IF NECESSARY 

1 Fli600.EQ»0) GO TO 414 
UG ^18 L^l tl^j 
liO 417 H^l ,1'i 

KHCS^(M)^KhCSZ(M)+GASP(L,Hi 
. 417 MiRoZ(L)=KhKSZ(L)*GASP(L,M) 
^Id i^Rf^:>«NR«S+KHKSZ(L) 
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DO '^16 M=l,I5 

C 

C THl:, lAI^GE matrix IS NECESSARY BELAUSE OF UCCUPAWCY KATE REQUIREMENTS 

rtwKH (Lf M, J ,K)=MWRH(L,M,J,K)+OASP(Lf M) 
^io oASHlLfM)~0 

lBLOKL r Jff K) :*LBLDt,(Lff Ji K )^IQLD(L) 
iULUlL)=0 

C 

C :50 rOU WON'T LOSE TRACK OF WHAT YOU HAVE ALREADY 

IF (iAV,RHD(L ,J f K) .GT.LKEQL(L , JfK) ) LREUL ( L i Ji K) = KWRdO( L , J, K) 

NCMH( J,Kr L)=NE*^RH( J,K,L )+MWKHD(Lf J,K) 

MWH.iC(LtJf K) =0 
413 lONTINUE 
^»J0 LCNFINUE 

iF(i I .£0.2) GO TO 430 

C 

C THIS CODING IS EXACTLY PARALLEL TU CODING STAKTINi AFTER 4 10 • BUT 

C NOT BY OISCtPLlNL AND LEVEL SU^JCE CLASSKUcJMS A^E INTcRCHANGfcAliLE 

CALL NLSf'I L( I BLD 1 13 t WRHS i NLREQi GASP ) 

iFl:>400.tO*0) GO TO 419 

DO 421 M = I,I^j 

DO 421 L=lfI5 

NHCSZ<M)=KHCSZ(M)+GASP (L iMl 

nHKSZ(L) ^KWSZ(L)+OASP(L iM) 
4^1 iMKMi-=NRMS*GASP(Lf M) 
4iy uO 420 L^I»15 

HftH:>LM(L)^URHSUH(L)+WRHS (L) 

NLBLC(LJ=NLBLD(L )+lbLD(L) 

iULO(L)=Q 

iF(,iRhS(L) •GT.NLREQ(L} ) NLREQ( L ) ^yJRHS ( L ) 

wUHi(L =0 

DC ^20 f' = lfll> 

wKHiUrt(L fM)=HRHZUM(L,M)+GASP(L,M) 

uASP(LfH)^0 
4ii0 CGNriNlui 
430 LUNriNUE 

C 

C TABLES 40If402 

iF(6^CQ.c'i.O) GO TO 300 

C 

C THE WEEKLY CLASS HOURS BY CLASS SIZE BEFORE THE HJUR SIZE FUNCTION 

C JK •SPILLING UP' ARE CALLED •BASE* IN THE TABLc 

riRIf E(N6, n6 )N,KHCSZ,NRHSjKHRSZ 

riKITh(N6,9'?9) kKViRH 
iuO LCi^flNUE 

C 

C rtHEWf "iHAT ENDS HOURLY LUDP, NUH TO PICK UP THE LOOSE ENDS AND MAKE 

C UP THE RtPORTS. THIS COULD PRUdABLY HAVE oEEN DONE BY SUBROUTINE 

C BUT THE* PROGRAMHLR HAD A PING-PUNG MATCH. 

DO :>C2 J=Ii33 

DO i02 K'^Lii 

DO :>02 L='IfI5 

AMSH( J <K »I n^WRHSAVlKf Lt Ji'^ACSlL) ^AWSH( JffKff i I ) 
•j02 uONF INUE 

1F(I I .EC.2 ) CALL PEEPl {N0i(5) iNCAMPfNYEAKiNITi ,NDATE , AI.SH, 2f 

1 ONAME) 
uALl KCLEA(^( UriliTQT,4i>) 

HHirECN6f 9131 NOU II*'4);NCAMPiNYtARiNPRa( i 1) iNl T , NRJNi NDATE, AFT 



00 tiOO J=«l |33 
UQ ^01 K«l,3 

CALL PlfEP^(JfWr<HSAV(i,itJ) lACSf lit IfONAMGJ 
t>JO cGN( INUC 

uALu PLep2(3^,h£fWTCTfACS tl 1 tOtONAHe) 
lF(i l.EQ.Z) CU TO 520 

hRlTctNat -n? ).NOi (II ) t NC AMP , NY t AR , NBGUM ( 1 I ) t i>i 1 T t N < JNt NOAT E 
rtRUti(N6f T10)DNAMHUf34} fONAME( 2f34 ) , NZE«0 .NLliLlJ 
WRirElN6f 919) NOi{ li ) f NCAMP,NYEAK,NIT,IMU*N .NOATd i NROOMi 1 I) 
CALL FCL(iABC0ATA,3) 
LALl ,KKi:ttR(NDATA,3) 
iWLA:>rT = J 

X=G* 

t;N=Q 

FPWSH^O. 

hPS^a. 

lFmLRtQ(L». tQ.O» GO TO 512 
NOAf AC IMNCATAii )+NLf<tQ( L) 

1 TEHpuNLcLCC L) *NLF^NS{ L ) 
hN=iTciMP*NSTACL) 

lTEi-»P-I feMP*FFPSNL(L) 
NWAoF ( if D^NWASF ( 1 1 ii +1TEHP 
WLA^F (LJ =Nl ASF(L) *-ITEHP 
NLAbFT^NLASF ri-NLASF(L) 
FPS^NLASFl L)/FLQAT(,NN) 
NDATA (2) ^»NCATA(2 )+NiN 

Y^Y/NLREq(L) 
uo :>ii f = i ,11^ 
iil X=X4-AC$IM>*WRHZUHIL,H) 

IF(X,GT. .COOOl) FPWi:H=NLASF(L*/X 

OATA (3) ^aATA(3)>X 

XaX/NN 

IFiY.GT. •0001) 2=X/Y 
oO TO 5li 
512 iF(WRhSUM( D.EQ.O) GO TO 510 

5iJ riKirE(N6,920 ) t NLf^EQt L ) i NN, WRHiUH < L) , Y , Z , X, FPS , FP HSH 

N0ATA(3):=NCATA( 3)+WKHSUM(L) 
51x> CONTINUE 

X==Nl)ATA( U 

DATA(i)=NaATA{3)/X 

DATA (3) -OATA( 3 )/N0ATA(2) 

UATA (2)=QATA(3)/0ATA(1 J 

FPS^FLGAfC NLASrT)/NDATAl2) 

FPhbH=FPS/CATA(3) 

iNASFdf I) »=NLASFT 

WRHEiN0f92l> NOATAfUATA.FPSfFPHSH 

oQ TO 200 

f>4L0 fciKlf l;(N6,9 17)N0l(l I) i NCAMP .NY t AK , NKUUMt 1 i) i i T , NrtUN , NOAT C 
wRlfElNti »922) NZErtU,N2eK0 
UQ :>21 J = lf33 
oa :j21 K = l,3 

i-il WRIT E(N6,9ia)0NAHEi iiJ JfDNAMEC2iJ) f K i ( LOIDL ( 1 1 J ) ,t^l,15» 
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>iRir E (N6 t919)N0l( I i ) , NCAMP ,N YEAK ,N1 T, NRUN , NJ ATt iNRU JM( 1 I ) 

LAU KCLnAP(NNHCL0f45) 

uALL KCLd/vK( Mi^^:0t^5) 

V.ALL KCLEAR(NCN-<H,45) 

cAU KKCLt:R(L»^eQf43) 

UU j JC J^L ,33 

UO !iiO K=»lf3 

CALt. FCL£AR(0ATA,3) 
CALL KKCLER(NCATAf 3) 

Y=C. . . 

^ = C, 

IF (LHCQL (L ,J,K) •EQ.O) GO TO 533 
WuAf/S (1 )=NCATA(l)+LKtQL(Lf J,K) 
1 TEHP--LOLJL( L, J,K)*LFNSlLf JiK) 
NN= lTt^P + NSrAL(LfJfKi 

i TlHP-I TEMP<^HFPSL( Lt J|K) 
NWASF (2, J) -NWASF(2f J) + IT[:MP 
LASF (Lf J fK ) = LAS»MLf J,K] + 1 TEMP 
lNA^F(2f J) =^INASF(2f J) + LASF(L, J,K) 
LReu(K,L) = LKtU (K,L )+LASF (L f J ,K) 
NFT~NFf+LASr (Lf J,K) 
NUAT A(2)^NCATA(2)+NN 
Y^NuV;HH(J,K,L) 
Y=Y/LRECL(L,J,K) 
UU :>i2 M = i,L 
t>i2 X^Xi-KWRHCL fM, J,K)*AC$(M) 
JATA (3)=0ATA13)^X 
X=»X/NN 

r4NHJLU(Kf L)aNNHQLO(K,L)+NN 
HHEJ(K,L)=NREQ(K,L)+LREQL(L, JfK) 
FPS=LASF(L f J,K)/FLUAT(NN) 
IFiA.GT. •0001) FPWSH==FPS/X 
1F(Y,EQ.0. ) GO TO 533 
Z-X/ Y 

533 r<UN.<H (K»L) =NONRH(K,L)+NtWRH( JfKfL) 
NO AT A (3)=NCATA( 3 ) ^ NE WRH ( J , K, L) 
IF (LREQL(L ,J,K) •NC •0,UR, Nh WRH ( J , K , L ) .NE • 0 ) 
IWRITt (N6,92J) ONAMf:! I,J ) ,0NAHL(2, JJ iK,L, 

2 LKEUL (L, JiK) » , NEl^ RH IJ , K ♦ L ) t Y t /.i X t FPS , FP'rtSH 

531 uQNT INUE 

iF(j>iOATA(l).EQ,Q) GO TO 530 
UATA (1) =NJATA(3i/X 
UATA13)=DATA13)/NDATA(2) 
jATA(2)=0. 

FPS^FLUAn NFT)/N0ATA(2 ) 
FPhaH^Oo 

iF(JATA(3),GT, ,0001) FPWSH=FPS/OATA( 3J 

IFUATAd) •EQ.O,) GO TO 534 

JATA(2) =UATA(3)/0ATA( I) 
53^ WHIT E (N6,924)CNAME( I, J J'f ON AME ( 2 , J ) t K, NOA TA , J A f A ,FP S ,F PWSH 
530 uUNflNUc 
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5^0 rtH [ rt (N!it9ifi> 5K , INASF( 1 tK» t INASF (ZtK) r NW ASf- ( I »K) , ASP ( 2 , K I 
wftir ECNOfSn INOKl I >tNCAI'lPiiNYeAK,NiT,Nr<LIN,NJArtfNr''^3M( I I , 
00 :;50 K=l ,3 
CALL FCLEAK(0ATA,3) 
oAt.L KKCLEi^i Ni;ArA,3J 

00 tfSl L = lfl5 

if (artGQiK, L) •ECO) GO TU 551 
NFT=NHT+L«EQ<K ,L) 

( 1)^NCATA{ l)+HREa«K»L) 
tlD AT A ( 2 } =N E ATA < 2 ) +NNHOLO ( K , L J 
NOAf A (3 ).=NCATA ( ^} -H^ONRHi Kf 
r--^NuNKH(K, LI 

ir=Y/M^ec(K ,L) 

A=NJNRH(K, LJ vACSlL) 

DATA(3)^iJArA<3J*X 

A=X/NNHULO(K,L) 

FPS=FLOATC LHEQ(!<,L) )/NNHQLO(K,L) 

IFU^GT* •OOOi) FPWSH=FPS/X 
iF(ir.EQ.O. ) GO TO 552 
^ = X/Y 

5t><: wiUVE (N6,926JK fL,MREg(Kf L) ,NNHUL01K,L) iNCNRH CK , L ) # Yf Zi X, FPS, FP WSH 
5i>l uUNr I^UE 
A^NUATA(l) 

iFlX^EU.O) GO TO 550 
UATAl 1) =N0ATA(3)/X 
J/iTA(3):=GATA(3 )/NDAT.;(2) 
JATA(2)=Q. 
FPWbh=0. 

rPS^f LQAT( NFT)/NDATA(2 J 

iF (JArA(3) .GT. •OOOi) FPWSH=^FPS/OATA( 3 I 

iF (JATA{ I) •or. .0001) DATAC 2)^DATa( 3)/QATA(I) 
5l*0 I'.RITE (N6f<^27)K f NCATA,DATA,FPS,FPW^H 
ZOO tONflNUe 
100 CONTINUE 
'900 FCRHAr{2I9,3Il) 

901 FOr<.^AT(//i i0=«fl9t' 1X=«|19) 

902 F0Pi1/ST(^A^f,I2,I4) 
90J t-URtlATi 15F4.0) 
90^ FQR/^AT( 151 3) 

90:> F0MAT(T67 ,12, I U Tl , i;)F^.3 ) 

90o irQRHAr(FlO.O) 

9U/ FCHilATl iOF7*6) 

90ti F0MAT(T7Iti;i,TitI5F4.3) 

909 FGRMATi I^X,3 112) 

910 F0f^i-iAT(l5i-<) 

9 11 rCPJAT CI2,To9 ,2 ( 2 1 T13 , 31 12 ) 
9t2 FURj-!AT(T69,2I2f Tl, 1514) 

913 rCKMAf (• iFAfI TABLE 2",A2,' FOK •fA4»» IN YEAk 'tUfM", 

lA3f«LAU tJEEKLY KOaM HOU^S tiY DISCIPLINE •/ < 11 LK A T I CN ' i 12 1 • , RUN 
2A4i't • f 2A4, OX'FQR EACH CLASS SIZE ANU COURSE LcVGLV/ 
JuXTiiESL- WHEKLY ROOM HUUH TOTALS AHE 1HE RESULT OF A GAUSSIAN'/ 
tuX'KANOOMUAriGN OF TKE CLASS SIZE U 1 b TK I UU f 1 UN • ,2 A3 , • DISfRlbUTIN 

VoX' CLASSES rO MUUKS ANO THEN 10 RUCMS ACCUHulNG TO OSO*// 
tj^7X« liEEKLY KCUM HOURS IN EACH CUUi<SE LEVEL • / 'OU I SC IPL liNE CLASS SI 
/^.E'f 7X« I'l IOX'2' ilOX«3»,^^X'TaTAL*//) 

9it iURiiAT { • IF/M TAULE 4 42 , /2GX • WEEKLY KUJH HJUi^S I cACH liOUi^»/27X 
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i'FOK CLASS • .A4f •S*//2X»H0URSi2X» i 2 J 4 5 

^ o 7 a 9 10 ii 12 13 14 16 TOTAL 

3 M 

yiu FQH/lAT(lM0,"eOA0, I SPUN OUT CF LUOP 41^'} 

9io FOHMAK 1H0,I3.3X'CLASS S UE; • , i6 lo ,/7X, •t^bC^ S 1 Zd : • i iX , 1 51 6) 
91/ FURrtAT(>lFAM TAJLE !>«iA2f« FOR »fA4f» IN YEAk •,l4,«:*i 

1* i^CV^LY CCf''PLeTtO CLASS SA4, •S-/» I TtKAT • i I 2f • HJN »fA4,»i 

i^A4//l3X-SU t: 1 2 3», 

d • 4 5 6 7 8 9 10 1 1 12 13 

3 14 15" /IX'UlSClPHNe LEVEL*/} 
91ti FORMAT (3X2A4f4Xi Lt^X15I6J 

')IJ FOHrtATC'lFAM TaBLF^ 3»fA2f' FOR 'fA4f« IN YIZAK •il4i': WEEKLY ROOM 
I ftCOM HOURS AND UflLlZATIGN HATf:SW* irEiUfiJN 'tI2i* RUNSA4, 
2» » •,2A4,YXMN CL AS f A4 , * S • / // 9X" THt ^jfc NUt^^itKb JF rtUh.vLY KOOM HUU 
iKS ARE AGCREGATICNS aF»/9X'THE FINAL NUMticRi IN ALL HUUKS IN EACH 
tKUQH SUE' /9X*VlITH ALL OTHER DATA COKKESPUfiO INC-LY AsiGKEGAf hO • // / 
D^OX'KCGM NJ, TOTAL TOTAL WKH PER UNi^TJ v^SH PER ASF PE 

oK ASF PERV DISCIPLINE LEVEL SIZE RCP:iS STATIONS WRH R 
oUQM STN OCC STATION STATION WSHV) 

9^0 FCR^HAT ( 7X« ALL ALL* f2I7fiafI9»Fa.l,F9,3,3F^.2) 

9^i FGRrtAT { 7X« ALL ALL ALL" i 16, " *• , 17, * v>c« , f j , » .f^ 7, 1 , • , P3,3, f « • , 
lJ(Fd,2,»'>' )) 

^^d FCRiiAT(2I2» 

92 J FGRi4AT { 3X2 A4 , I 4»2 I 7 . 1 3 , 1 9 . FO • 1 , F9 • 3 1 3F9 • 

924 FGKi4AT i iX2A4, 1 4t4X, 'ALL* il6,«««,I7,"*",f8f*<^*iF7-i,«*',F8.3,«*', 
i3(f3,2,'*» )/) 

92!i FORHAT( 12|4I10» 

920 F0PiiAT<7X-ALL« ,I5,2l7tI8,I9,F8,l,F9,3,3F9,2J 

92/ FORMAT (7X» ALL* ,I5|4X»ALL" t l6,'»M7,"*Md,'**,F/«lf«*'fF8.3,***t 

i3(Fd.2t )/) 
9rl<i FQPMAr(T70tI2iTl,3FL0.2l 
^29 FORMAT I 20A4) 
999 FOP.MATC 7X, •BASE' t7Xt 15U) 
tND 

SLtiRCLTINE HRS ZF N ( KKWRHt ACS, JWRH , KWRHDi IXt 1Y,FX) 

iNTcGER*2 KKWRH( L5> , JWRH ( 3 , 15 i 3 3 » t KHKHDI 33 ,3 , 1 5 ) 

uIMcNSlON ACS(15) 

C 

C THE FULLOWING LOOP THROUGH STATEMENT NUMBER 12 HAS The SOLE PURPOSE 



C Of ADJUSTING THE NUMBER AND SUE OF WlE^^LY ROJM HOURS UNOcR 

C ThE INFLUENCE OF THE CLASS SUE FUNCTIGiJ. WEEKLV STUGENT HOURS 

C BY HOUR, DISCIPLINE AND LEVEL ARE KEPT APPROXIMATELY CONSTANT. 

C TO DO THIS THE FOLLOWING STEPS ARE PCKFOKMEDi 

C U THE CLASS SIZE FUNC TI 0:4I FX) IS FOJNL}. i MAIN ROUTINE) 

C 2. FOR EACH CLASS SUE A NEW ACTUAL CLASS SX2t IS FGUNO BY 

C IPLYING THE OLD NOMINAL CLAi>S SUE liY FX. 

C 3. THE STUDENT POPULATION IN THAT HGUR AND N]MINAL CLASS 

C SUE IS REOISTKIbUTtU TU THE TWO CLASS S12CS(0F THE RANGES 

C AVAILABLE) ON EITHER SIOE OF THE ACTU^^^L CLASS SIZE THAT 

C HAS BEEN FOUND* THIS IS DCNc IN SUCH A JAV AS TO TRY TO 

C KEEP THAI ACTUAL CLASS SUE ANO TJ KEEP IHE STUDENT 

C POPULATION APPROXIMATL-LY CONSTANT. 

C 4. THE NEW NUMBER OF WEEKLY ROOM HOURS BY MQ'JK AND CLASS 

C SIZE MAY UE BIGGER OR SMALLER THA.J THE ULO NUMBER. JN 

C EITHER CASE THE NEW WEEKLY RCCM HOUS^.S A.< ti OISTRIBUYEO 

C ( PROPGriTlONALLY TO THE OLDiACROSS JISCIPlIWE AND LEVEL. 

C 5. IN GENERAL IT SHOULD rtE POINTED OUT THAT THE NUMBER OF 

C iiPECIFiC CLASS SUES LIMITS THE Q.JALITY OF THE ROUTINE, 

C 
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if (KKWKHCLLi -Eg. OJ GO TO iZ 

if lAtiS ifX-i.) .LT, .01) GU TO 512 

if ifX *GT, !• ) GO TO oi2 

C 

C THi5 li UHcN CLASS SIZE WILL OECftHASE 

C LLL WIlL point to HiGHEX OF 2 NiZW CLASS SIZES 

C KKWKH IS NUHdER OF ROOM HOURS OKIGINALLY ASSIG.NEO bV SIZE 

LLL=LL 

bO TO 1103 

C 

C bJilP OUV.N UNTIL COKRECr PAlft Of^ CL .S SIZES IS fOUNO 
iiui. IF IFX .GE. ALS(LLL-i) /ACS(LL) J GJ TO IIQZ 

LLL=LLL-i 
iiiJ3 IF (LLL .NE, I) GO TO 1101 

C 

C IF WE AKE Af 8UTTCM OF LIST ASSIGN EVERYONE TO RIGHT NUMbER OF 
C brtALLEc^T CLASSES 

.>iA = l KKWRHl LL)«ACS(LL) )/ACS( i )*«5 

LLL=2 

GO rO 1106 

C 

C FKUH 1102 TO JUSi BEFORE 512 IS ROUTINE TO FINO ACTUAL NUMBER OF 
C ULAoSEc> OF TWO APPROPRIATE SIZES HAVING PINPOINTED THOSE SIZES (LLL) 
1102 A3=»FX«ACS( LL ) 

P^^AOS (LL)*KKWRH ILL) + 01 FF 

IF I H ,GT, 0. ) GO TU 1115 

UlFr^P 

UO TC 11 

111^ V Y^P« CA3-ACSCLIL-1) )/i AS^tACSfLLD-ACSiLLL-'l))) 
Nb=Y Y-t.S 

NA«( P-A3*YY)/A3*,5 

C 

C HAVING FCUNO THE NUMBERS OF CLASSES IN EACH CLASS SIZE WE TRY TO SEE 
C IF THE LOWER CAN BE CHANGED TO A BETTER FIT (BECAUSE QF INTEGER 
C WOMbER PROdLEM:^) 

NCUM^'NA^ACSCLLL-D+Na^ACSCLLLi + lNA+NB MA3 

XX=( 2.*P-N0OM) /(2.<'ACS(LLL-1) ) 

NA^NA + XX + .5 

IF (NA .LT. Q) NAaO 

GO rO 1106 

C 

C FX iS TCC S^ALL TC BOTHER TO CHANGE 

512 UO 51J K^'l f3 
00 513 J=^lt33 

JWRH(K» 16-L» J) -KWRHOt JfK,L) 

513 KWRHCi JtKf Li=0 

GG rc 12 

C 

C THIS IS ViHEN CLASS SIZE IS TO INCREASE 
612 LLL:=>LL + 1 

oJ TO 1105 

C 

C 3UMP UP UNTIL CORRECT PAIR OF CLASS SIZEb IS FOUND 
llo^ IF (FX *LE. ACS(LLL)/ACS(LLn GO 10 1102 
LLL^LLL+1 
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C iF mE are at top of list assign EVr:RyG*s£ TO KIoHF NU.IkJiER OF LARGEST 
/ C oLAbSES 

1103 il' ILLL .NE, 16> GO TO 1104 

t^A^i KKWHHUL}*Aci(l.t ) ) /ACS ( 1 5> - 5 

C 

C TrtU NOOLi WHICH SUURaUTlwE MUST GO THROUGH 

C IF is-tAoSES TO BE ASSIGNLDt CLASS SIZE PAIR IS FJUNO ANJ NUMBER np 
C IN EACii HAS BEEN CALCULArED 

C AOUM A-^U OIFF CALCULATION iS TO KEEP RUNNING TOTAL OF r^UUNUlNG ERROR 
C UF hSH and ADJUST FOR IT IN NEXT LOWEST ROOM SUt 
llUu AULM=NA*ACS{ LLL-1) 

AFiNC .NE# 0) ADUH=ADUH4NB*ACSC LLL) 

OiFl-=iACS( LL) *KKWHH< LL > ^Oli'H^AOUM 

IF 4NA+Nb .EU. 0) GU TO 11 
C XX WILL 6E THE DISTANCE BETWEEN UKIGINAL WRH THAT WILL bE CHIJSFN 
C FUR TH£ TOTAL UF THE TWO CLASSES. 
C WILL 86 T^<E DISTANCE BETWEEN 

C iHt CHJSEN WRH THAT WILL BE CHOSEN FUR THE HIGHER OF THE TWO 
C CLA^SEo. if MUST bE 1 OR GREATER. 

C IF i4B'=J THAT MEANS ALL v^KH Art TO bE ASSiGNEO TU LOWER CLASS SIZE* 
C PulfifL IS POl^TER TO THE NEXT ORDINAL HRH TU BE CHOSEN. 
C f'OiNTH IS PUINTER TO THE NEXT ORDINAL CHOSEN WRH TO dE USED AT 
C HloiiER CLASS SIZE 

V Y=NAf ^d 

POIWTHcO. 

AX^KKWRh(LL)/YY 

IFINB .EQ. 0) GO TO 1107 

YY=YY/N15 

C KA.^jGM NUMBER CALLS ARE USED TO MAKE STARTING WITH FIRST WRH A 
C KANOCH EVENT 

v,ALl RANOUdXf lYf DUMH) 

1X=1 Y 

PQlhlTH = 0UMH«YY + .5 
110/ lALl RANOU(IX»lY,OUMM) 
1X = 1Y 

PUINTL=UUMM*XX+.5 

C NTUTE ANO NTOTE^ ARE RUNNING TOTALS GF ORIGINAL WRH AND CHOSEN WRH 
C SO FAR CCNSiUEREO. 
UlJ NfUTE^O 

iucrE2-o 

C NPr ANU NPT2 /RE INTEGER ROUNDED VALUES OF PUINVL ANO PUINTH 
C FUR COMPARISON WITH APt^kOPRIATE NTOTbS* 

rJPT-?C INTI.+. !> 

NPT2=PCINrH+.5 

bU ilOa J=U33 

UO IlOb K^l,3 

if I K\tRMQ( JfKtL ) -EQ. 0) GO TU 1108 
NTUrE=NTOTE+KWRHD( J.K, L ) 
NWRr(0( Ji^f L) =0 

111^ IF iNPT .;;t. ntute) go to iloa 

IF INB .£0. 0) GO TO 1112 
NraTE2=NrOTE2+l 

iF (NPT7. .GT. NTOrEZ) GO TO U12 
L^t^LL. 

PCIliV =P0' NTH4'YY 
WPl^=PaiNTH^.5 
uC TC 1113 
1112 c-^^LLL-l 
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iliJ PCINTL=P0INTL-XX 

jWfid IK,L4, J) ==JW»<H(KiL^, J ) +1 
rC 11 H 

lloa ccNf iNue 
GO rc 12 

ii OQ 10 J»lf33 

UU 10 K=^l,3 
10 NWf<r<0( J|K,L)=a 
I*: CGNf INUE 

K6TJRN 

jUti^CUTlNf: KCLe>^R(K,N) 

1 K(L)^0 
KtTUKN 
tNO 

bUBKCUTINE rtANOUdX. I> ,YFL) 
C i>fc6 SClcflFlC SUBROUTINE PACKAGE - FLAT 0ISTK16JTIGN RANDOM NUMBERS 

4 Y^l 

irCiY) 3f6,6 
3 i Y=l Y+2147A83647+1 
o YFL-^IY 

YFL=YFL«,^6i)6613E-9 

RETURN 

ifND 

StBHCLTINfc GAUSSdXf S.AMfV) 
C SEii SClENflPIC SUliKOUTINE PACKAGE - NORMAL 0 IS fi^ I BUT I ON RANDOM NUMbERS 

00 jC 1^1, 12 
CALL RANUUdX, IY,Y) 
IX^IY 
5j A==A<-Y 

V^(A-6, ) 4S+AM 

RETURN 
cNO 

oUORCUTlNE NLSPILH8LD,SS.NRf NLNR.MWRH) 

C 

C fHlS SUaRQUTiNE ASSIGNS CLASSES ru CLaSSROCMS* IT STARTS WITH THE 

C LARGEST CLASSES, UHLN THERE ARii NOT tNUUJH KUJMS THAN A LARGER 

C ROOM IS SOUJHT, IF THERE AKE NUT ENOUGH LARiER RO JMS THAN A ROOM 

C IS BUILT, 

C 

C .^R IS THE NUMdGR OP CLASSRUOMS NEEDED OF A GIVEN S(Z£ 

C i^iLNR IS THE NUMBER OF AVAILABLE RUCMS CP A .HVLN >IZE 

C 13LJ IS THfc NUMBER OF ROOMS TO OE bUILT OF A GIVEN SIZE 

C HLhA IS THE NUMdCR OF EXCESS ROOMS OF A GlViiN SUE 

C HWRH IS AN ARRAY OF CLASS SIZE BY ROOM SUE SINCE CLASSES MAY BE 

C ASSIGNED A ROOM SIZE NLLD LARGER THAN THE ACTUAL CLASS l2r.rORE 

C cNTRY UTO THIS ROUTINE (THROUGH USB), 

C K IS THE RCQM SIZE 

INTtGER SS 

lNTuGER*2 IBL0(l),NLEX(15),NR(I3l.NLNR(i)iNWRH( 15, 13> 

CALL KCLEArt (NLtX,lt>) 

K^^SS 

C 

C FIR:>T FOR EACH SIZE FIND THE EXCESS U>s SHORTAGE OF ROOMS 

iF(iNR(K) .GT, NLNR(K») I bLD( K ) ' NR( K )--NLNRl K ) 
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IF IMNR(K) •GT, N^xtK)! NLEXIKI - NLNR IK HNi< ( K) 
to ^ = K-1 

if IK .EO. 0) retuk;h 
JlTfO 

iH-(.^K(K) .GT^ NLNR(K)) I BLO ( K ) =NR ( K )-NLNR i K ) 
IF (NLNK(K) eGT. NR(K)) NLEX(K) » NLNRl K) -Nr< CK ) 

IF ThEKF IS M SHORTAGE TRY THE NEXT SMALLER SIZE 
bO if (leLO(K) .EQ* 0) GO TO 10 

If THERE IS A SHORTAGE SEE IF A LARGER RCUH IS AVAlCAOLE 
if (NLe)^lJJ .GT, 0) GO TQ 30 

iF(J-SS) 40,^OtlO 

IF A LARGER ROGM IS AVAILABLE ASSIGN If AND aOJUST ALL VARIABLES 
30 IBLUIK) ^ lOLDiK)-! 
wLEXt J)=NLEX(J} -I 
NKCJ)=NR( J) + 1 . 
tiRU) * NR(K) - t 

MNU THE Lyi^AGEST GLASS SIZE IN THAT ROOM SIZE AND REASSIGN IT 
KK IS THE CLASS SUE 

7u IF (f<kRH(K,KK} .EQ, 0) GO TO 60 
rtWRH( J iKK) =MWRH I J, KK) 4-1 
rtWKrt (K,KK) -M>JaH(K,KK)-l 
\jQ [Q 50 

IF(KK .NE, 0) GO TO 70 
HRirE(6y 100) 

lUO FUHi4AT(» IMPOSSIBLE tRROR- I 
kETURN 
tND 

SliUKCUriNE FCLEAR (A,NJ 
uIHtNSIGN AC 1) 

OU i 1^=1 tii 
I A(L}^0. 

REfORN 

bUliKGUT INti KKCL£R(N»;.) 
Uiry^NSIQN Hii) 
OQ i K=»l ,L 
i i^j(Ki=0 
KETJRiN 
tiNO 

bUBRGOTINt ERRaR(N^RK,NI 
iSiERK-NEK^^-l 
wRUE (6»]LU0) NfNEKR 
lOU hOMAT(»OEHRGR NUMbEK 'tI2t' FUK A TOTAL UF SI^* • ERRORS, PRGGRA 
IN HiLL STOP AFTliR READING THIS GROUP UF INP JT CARDS, M 
KEfORN 
tNC 

SUaRCLiTlNE PEEPl(NO,NCAMP,NYtAK,NI TiNRUN,riOATE t WSH,HC» ONAME) 
UINL^SIQ^i CATAZ(3) tV;SM(33f 3i2) tOAfAii) ,N0< U ,NDATL(Z) 
Ul^cNSICN SUarOT , ONAHcC^rJ-^ ) 
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V*RirE(N6 f I 19) Na{ I) fNCAMPiNYEAR ,Nl TfNRUNfNCATE 
AFCi'lC.EQ.^ J WRITE (N6rl203 
1F(A1C.GC.2 ) GG TO 10 
V^KUt (N6f 1211 
lU hKir E(N6r 130) 

LALu FCLeAf<(0ATA2f 3) 
CALL FCLEARt SU8rOT,8) 
UG ^60 J-lf33 
CALL FCL£AR(UATA,3i 
X = 0. 

OATA (K) ^OATA(J^) + WSHUf K. ! > 
6UBrGT(K f I ) = SUGTUT{ 1 )+WSH( J,K, I J 

:jUb J 01(4.1 )=sui3raTi'»,i 3+*r^SH(j,Kf n 

2o2 Y = Y*-;>SHl>Jf K,I) 

iF(r .N£. 0^) WrtITC(N6f 122) 
lONA/iEd ,J) ,DNA«i:(2, J) f 1 t (WSH( J,K , 1) ,K=if 3),^ 

26i x=;<+Y 

IF tX .N6. 0.} WRnECN6, i23)UNAHE( 1, J) ,0NAMc(2, JJ , OATAfX 
^.=Z4-X 

UATA2(K) ^ QATA2{K)+0ATA(K) 
ZbJ CONTINUE 

hRirE(N6,124J SUBT0T,DATA2ti 

RETURN 

Id^ FQRHAT(3X2A4» 9X, 1 1 . 17X, F9U t 3XiF9 .it 3X,F9. I, 5X,r-9. I) 

i^^ ^aR-^A^ ( 3X^A4i27X,F9.I, t2X,F9. i2XfF9.li '^^^ ,^X,F9, J., 

liL^ hCRi-lATtiX, »ALL* t i2Xf * 1 • t i5X» 3 ( ^X, F9 . 1 • • * • > i ^Xi FS« I , • * » / 

1 iiX, 'ALL' ,12X, , l5X,3(2Xt 1-9.1 f f4XiF9,W»*V 

2 5Xi 'ALL S 29X, 3{ 1X,F9. If ) ,3X,F9*lf ) 

11^^ FQRi^ATC IFAM TAauE lSA2i« FOR SA^f' IN YEAR. S^^f «: V^EEKL 

lY STUDENT HOUi^S IN EACH DISCIPLINE*/' I TEf^ AT lOii* , I3f • OF P,UN 
2A4f> ON •f2A4,5Xt*0Y CLASS TYPES IN EACH CJUKSE LEVEL*/) 

JL2i FQRMAT(/6Xr* THESE* , 

3 »^UMi3EHS OF WEEKLY STUDENT HUURS AfU- THE AGGaEGAVEa PRUCUCFS OF*/ 
^oXf'THE PKGJECTED ENROLLMENTS OF EACH STUOENT LEVEL AND I^AJOR* /6X 
'i>*ANJ THE iNUUCeO COURSE LUAO MATRIX (NCNLA6 LAJ HiH PER STUDENT* 
o/ ) 

12J ^CRMAT(/6Xf* THESE ARE THE FINAL AGGREGATEJ NUMbrIRi OF WElnLY*, 

1* STUOENT HCU.^S*/6X,*bENERATE0 OY FAMSIX.*/) 
iJO ^0H.-1AT(// 

b • DISCIPLINE CLASS TYPE CGUkSE LEVEL: I • t 1 IX, * 2 * f ilX , • 3* , 
rJX,* TOTAL* /) 

tNO 

iVJBRGUTINE PEEP2 II , JWRHfACS , U # KSWTCH , UNAME ) 
U1M.:NSI0N ACS I li ) , NOATA { 4i tAVG(4) t CNAME < 2 , 34 ) 
lNTcG£R«2 JWKH<3,15) 
Nt)=2 
U6 ^ 6 

CALL fCL£AR( AVGt4J 
uALl KKCLER(N0ATA,4) 
JO i L^liltj 

H^O 

JC 2 K=>1,3 
M^Mv Jk^RH(K ,L) 

^^0/\r AlK>«NCATAiK>+JWRH(KiL» 



^ AVGl K ) = AVG (K)*AC5(L )*JHP.HiK,L) 
ityNAHH iifU tDNANE (2i i J iLi ( JW/<H(K,Li ,K«I»3) i f< 
NOAFA (^) ^NCATA(n+N0ATA{2)>NaATA (J J 

IF iriOATA(4) .eg* o) return 

hKlTE (N61IOI) NDATA 

iF(KSViTCH.EQ.l) WRITc(N5,103) H. I , (NCATA(K ) I, 3) 
AVOl ^)=AVG (L) *AVli(2)+AVC(3i 

lFii>»DATAU>.EU.O) NDATAlK)=^l 
i AVGlK)^AVG(K)/NOATA<K) 
WlUri: (N6,I02 ) AVG 
HtTURN 

100 fORHAJ (3X2A^t 7X t 1 Zt 8Xf 14 1 7X i { 7X , 1 4 , 1 1 j ) 

101 eORt-lAT •TOTAL WRH" t4X 1 1 4 t » * • 1 6X» 14 » ' « • t (>X , I • • , 6X * v • * M 
FURMATl l!3X,«MEAN SU t • ,4X ,F6* 1 1 3Xf F6* i , 5X t FO. 1 1 5X , Fd. 1/ ) 

i.U3 FCR.^1AT1212/3I4) 
tND 



72 



VARICOS - PROGRAM LISTING 
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r 
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r< 1 
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bK 
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lUii' iC <>M 
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i\ 1 
U JL 
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V* Mr i-» 3 I 


r 
\# 


I A k> ZiT 










r 
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L/r J r 1 


r 




r 




r 




c 




C 


FAO 


c 


FaCOST 


c 


FaCK 


c 


FACIAL 


c 




c 


INiUM 


c 




c 




c 


WJAT E 


c 




c 


NX r 


c 




L 




c 


NVcAK 


c 




c 


Urt:>uM 


c 




c 


UY 


c 




c 





YEAU 



ACAJtMIC STAFF 



S^'ACc COiT TOTALS 



4) 



CUNSTRUCTCM PERIOD 



ULCTIUtNAKr' OF VARiAt3Lt NANhS 
The Wur<L 'ALL • 
TA SHA'Kt GF WRH WORKLOAD 
SHARE OF WRH WORKLOAD OF 
TOTAL VARIABLE COSTS PER 
BOND LIFE 

aCiNU INTEREST RATE 
FACULTY COST SUUTOfALS* ALSO 
ViCH/FTl: FOR TA«S 

ViCH/FTt FOR MON- TA ACADEMIC SfAfF 

YEARLY COSTS OF CAPITAL CC^iST:<OCT 10 ^ FROM YEAR i START 
CAPITAL CGSf FACTOR: lNUEX=IOa 
DISrRIi;uriON Of CAPIIAL EXPt:^J^TU.^l5S OVER 
SUPERVISION OF TA" S; ACAUEMlL STAFF/fA 
DISCIPLINE NAMES 
$/ASF: overall AVERAGE 

FACULTY OVER^RlDt FACTO!': RtGULAR FACULTY/ STUDENT 

$/ASF FOR M&O PLA.WT 

ENROLLMENTS 

EQUAL PAYMENT FACTOR 

FTE FACULTY ( TA 1 1 HKEG. , NON^-TENURLO , T E.>JURED ( i--^ ) 
SALARY COST'S 

ACADEMIC STAFF LESS TA'S 

AVERAGE FACULTY SALARIES BY -'UNK" 

ASF inVtNTOKY:NLiN-LArifLAd rRESEAi<CH, JFr ICEi SUPPORTING 
YEARLY TOTALS: INVENTORY ASF 
LCOP C'viNUOL FOR FINAL TAliLc PKl NT-OUT 
CAMPUS i\AM£ 
CALEhOAK DATE OF RUN 
TEMPORARY VARIABLE: NON-CLASS 
ITE:?.ATIUN NUMBER (N-TH YEAR OF 
NUMBER OF RUN 

NUMBER OF YEARS FOR THIS RUN 
YEAR OF DATA OISPLAYEU (INPUT 
NEW ASF ;NON-LAB,LA[i,RLStAaCH, OFFICE t SUPPORT iUG 

YEARLY TOTALS :NlwlY FINISHEU ASF 

TEMPORARY VAR TAULl : I NVE NT(jR Y CF UUC. ASF TOTAL 
TEMPORARY variable: Ntr^ DISCIPHUC: ASF lUTAL 
OTHER OPERATING COSTS YEAKLif TOTAL ( S AL AR I ES ^SUPPOR [ 
OFF ICL ASF/FACULTY 



REQUIREMENTS 



AS FI.^ST YEAR OF STUDY) 
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SUPPORT RATE: $/NGN-TA ACADI:MIC STAFF 
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SUPPORT ASF/RfcSEARCH A;4U UFi-ICb 1 U 1 AL 
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SUPPORT ASF/LAB TOTAL 
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TEMPORARY VARIABLE: iTE TA'S 
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TEMPURAKY VARIABLE: bUM OF FACULTY AND TA SALARY COSTS 
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TEMPORARY VAKI ABLE : TOTAL FACULTY iiUPPORT COSTS 



iNTuGER BL 

INTcGEH V^'RH{2,J,33) 

lNTcGER-*2 NARK(14,4> 

JI^<cNS10^i ALL(2),Ai(2,3f33),A2( 2f J f 3 J ) , BUCKSMi o ) , -UC 2f 3f 33 ) ♦ 
i 1)2 (2 1 3, 33) tCAPCSTd^! tCAPF AC 15*33) ,CSTUST (" ) ,C2(2.33) , 
2UNA{'1L(2f33 )t02(4t33) f tNR(4,33) ,FACl't,3J) ,FACUSr (5) tFACR( 33) , 
3*-ACSAlK) , l*NASFi5,3J) , 1 NSUiM { 6 J , NOA TE I 2 ) jNE b US 1 i ), N?^A SF ( 5, 33) • 
tNWSJ^(6) fOOPS ( 10) t UK 2 ,33) t PLOP { 10) .PRJCSTto) #02(2 ,33) , 
i>KK01ST(3t33J ,RKSUM(b) ,R1 ( 2 . 33 ) , R2 1 2 , 33 ) i SUPRA r(J3) . 5 H 33 ) , S2 ( 33 ) , 
oTA(2) 

OATA ALU L)f ALL(2]/'ALL 'i* •/ 

iN2^L 
10UT=6 

CALu FCLEARiCAPCST, I't) 
CALu hCLEARl iNASFtibb) 

REAiJ I INit^)0O)NCAMPf NKUN,NL)ATEtNUMYR^,NYEAR 

J.N=.JUWYRS*'4 

Rl:AU (IN2ty0i ) ONAME 

00 JlO K=Lf33 

R6i5J( IN2t90't) ( (Al(ItJiK),J^I,3),('=I,2) 
KCAU( IN2,90<t) ((A2(IiJ,K),J = l»3),I = i,2) 
KE^j(iN2,9C^f)neL(ItJiK) ,J=:I,3) fl«l,2) 

10 KOAu I IN2,00-^^) ( (iJ2< 1 fJi K) ,J^I,3) , 2) 
kEAJ(1N2,^^0^) FACSAL 

UO il K--.il 33 

11 KEAj IIN2 »90^) (C2{1 ,K) ,1 = 1 ,2 ) f {02tJ fK ) I J = i,'t ) 
00 12 K=ilt33 

1-i rGAu( IN2,90'^) SUPRAT(K), (rsK01STlJ,K],J-i,3) 
C tl4L) GF INPUTS FOR FACULTY C GMPUTAT 1 GiNS 
UO 20 K = l,33 

^0 RCAu (IN2 , ;'03) RU If K) ,R1<2,K) ,R2( l,K) , R2l2tK J iUi( I ,K) ,01 (2,K) , 
1 02U,K) ,02(2tK) tSIlK) »S2tK) 

UQ 21 33 
21 KEAO(IN2,90i>) (C APFAC (I i K ) , I--- 1 , ii I 
2<: KIiAUllM2i^Ga) K, lINASFd ,K) ,i«3i5) ,NUCHK. 

iFli>li:CHK,EQ.O) GO TO 22 

KEAJ( IN2,007) SiNFL,faR|iJL 

KEAi> (IN2i905) CSTOST , OAMNO , UP ;QF f 
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KS=HS+CSr05T(5-NJ 

c 
c 

UO i^JCO Nl T =1 tNUMYRS 

LALL NCLEARC NViASFf 165) 

LALL hCLuAR(HUH»198> 

^.ALi. FCLliAHC tNKf 13^J 
IJi KfciAJ( INat906) K, (ENKd ,K] J 1 = 1 ,4),N0CHK 

IHd^CCHK.cQ.O) GO TO 101 
lU^ KEAJ( lMf902) i,K 

.CQ-O) GO TO 103 

K£AU ( 1N1,^?03) (VlRiU i i J i K ) ♦ 1 , 3 ) 

bO rc 10^ 

iOi CALc FCLlE AR1FAC,132> 
.AL*- FCLtAR< 6UCKSM,6) 
v-ALu FCL£ARIRKSUN,5) 

WRirt HGur f9 50)NY£AK,N0AT£,NRUNtNlT tNCAMP.NVEAK 

C 

UQ 1 10 K^l f33 

LALL FCL£ARtFACQST,5) 

FACl<lK) = Q. 

6UPCST=:0* 

UO 1.11 I=lf2 

rA( I >^o. 

DC 112 J^lfi 

IFiWHHli fJiKJ^EQ.O) GO TO 112 

I A(l ):»(A1( Ii JfK>»WRH(I ,J,K}/Ol(I,J,K) ) + TAtl) 

FACK (K^ =^FACR (K> + ( A2nf J,K)*WKHU» Jt K)/02( IfJfKU 
11^1 CONFlNue 

*^AC{ 1 tK)=r-AC(l,K) + TAin 

F ACR(K)=FACR (K) + 102(1 ,K)«TA( I ) ) 
lijL LUNFINUe 

UQ 113 Ul .4 
11 J FACr<«K)^-FACRU)*{02t If K)*ENR(I,K) ) 

X-^FACl 1»K)4FACR(K) 

lF(X«Lf.0.n GO TO 110 

ua 1 14 I ==2 t4 

U-^ FACi I ,K) =:FACKi K)^RKUIST{ I-1,K) 

5UPuS r = .su?MAr(K )*facr(ki 

DG 1 lli I-» ,4 

rACuSTC I )^FACSAL« I )<*FAC< 1,K) 
bUCK^W< I )^dUCKSMn )+FACOST(n 
KKSUMd )-RKSUM( I) «FAC(1,K] 
FACa:iT(f>)^FACaSr(iiUFACOSr(ii 
YV-SLFtSr + FACQSr 
bUCKSM( 6)=»eUCKSM(6) + SUPCST 

WRXrE( iaUT,9i31) (0NA:U:( IfK) ,I=l,2»i (FAC( ltK)f l^lf<r) tXtFACOSTr 

1 SLI'CSTtYY 

liu LC^r uue , 

un 120 V- lt4 

bUCrwSM( 5 )«0UCKSM(5)+0UCKSM( J) 
12J KKSUM(5J^f^KSU;'l(5URKSUM(JJ 
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LUPj (fi I T ) --cUC^:n{ 1>) +oUCKSH(6) 
v*Kll t; ( lO'jr r^J:>l) A'lUtK.\w'.'M,tJUCKSH,UlJPSlMT) 
C OH HACtLTY CGMPUTAftCNb - b T Art f :)PACt CALCUtAT I QNS 

1 KilAu ( IM f'UO J K r INASF (1 tK) , IN ASF [^,K) f NWASF( i.K) , AiWASFt 2,K ) 
U ALL FCLc AK^:)ljCKSMt6) 

CAlu ^(.Ltl AKl ir;SUHT6) 
L ALi NCLSAHJ NWSljM,6) 

WAIi L I ICUr f )o:J)NYtAR,NCATt iNRUNf.Nl T ,NCAMP,NYEAK 

C 

.>JY=J 

uALL FLLeA«(Prt jCST,ol 
CALL NCLtAR(NtEDSf3) 

lF(t-AI-;v(<J -L r.O. L.AKO.r:NK( 3,K) .LT.Q, I • AND, t :{K ( 'i t K) .LT.O. DGQ TO^IO 
Nttui,( U = U ( If K)*» ACK< K) ♦KKZfK) ^FACC i,K)+KJ(i,K) ^t-NKl 3. K) 

NEEUS U) =ai( 1 ,K)*FACKlK)+01i;^ ,K) ♦FACt I , ► ) +Q 2 ( L , K ) * CNR ( 3,K) 

uLf:Ot> C3J '-^S UKJ *(i^tLUS( 1) +HE£:DS( 2 ) ) +S2 I K . * I N A SF i 2 , j< | 
JC 2 10 1 ^1 fi 

lF(/itEOS( n .LT* lNASFU+2iK) ) GO TQ 210 
uWA:>r ( I*2iK)^r^eEUS( 1 )- INASF ( I*-2tK) 
XNAiF(i+2f K) =NCEOS( I) 

210 LCNr INUE 

UQ 211 1^1 ,5 

PrtJubH [ ) = N;^ASF (I,K)«CAPFAC( I ,K)^aR*OPSQFT*f^$ 

PrtJi-ST16' '^UCS^Idi)^^PKJCSTCI ) 

i^i^^-z^Z-HNASF ( I ,K) 

;4Y=*\Y + Ni-.ASF( 1 ,K) 

UWSJMd )=NWSUM(l)+NtJASF( I,K) 

iNSu^'i n^i ^suM( i)'^iNASP( r ,K) 

211 atCf<.SM( I ) =nUCKS«( 1 )+PRJCST( I ) 
IFINX-Eg.c; GQ TO 200 

riRr L (lOOr ,961) (UiNAMEd ,K) . 1 = 1 #21 , 
A ( INASF (l,K),l = l,S>) (NWASF( I .K) • I==l ,51 .NY, 

IPRJCST 

200 ugn;inue 

UQ 23C J=l,5 

WWSJ/^ 16) = NV>SUM(6J ♦NhSUM( J) 

1 NSOH (61=«[ NSUM(6)+lNbUK(01 

2 30 oUUSf.C 6 > = bUCKSMlo) +tiUCKSM{ J ) 

bO ^20 f 1=1,4 

XX=iiL0KSM(6)*EPF'XCSTOST( lU/RS 

iviu=;aT*i (-1 

JG I =riC» JN 

uAPv-sr ( I )==CAPCSTjn+xx 

220 LGNf IiNUrl 

riRiri:( IGUr #961) ALLi INSUM,NWSUM,bUCKSM 
PLLV (rn r ) = lNSUM(6)=<^U4HNU«f'r< 
l.\' ' INUE 

J = ^YEAR-3-^UMYR > 

WRIT t ( IGUr ,9 70)NOATr:,N)^UN,NLiMYRS 
uG iOOl 1=1,4 

CALL CCNVfl(NAPr<(l,I) , 1 1 1 C A PCS T ( 1) ) 

r,KIitClJur,S? J J .(NAkRUf 1) ,L=1, 1<*) , (NARRCLt 1) tL^l,i4) 
lOOt J^J+1 

JQ jlOC2 I "UrUHYRS 
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uALL cc^v^. (NARK ( , i 1 , plc;M i ) j 

I- ALL ^Cl^V'MNA;s^^( t ,ii,CAPCSr ( U-tU 
uALl- CCiNV.< (N.\kK (1,4) t i I? aiGSUM) 

HKH E: ( iG»jT , ; a) J,NAHR 
vuj i-crMi/^T('»A^,iaf i^) 

Slui i- OH.iA r ( 20A^) 

VU't HCPMAr(OFlO.O) 

9J/ i-CR.-lAK ^^5*UtI2) 
vud i^CK.-JAT ( 31 10, 12) 
910 f-Gi<.MAT ( 12,^1 10) 

9:^0 hQR.iAr( ' IV^^iCaS T'AOLE: A-« , I^V , /,a/ • CAT £ : »i^A4/* KUN NUMSliK: "lA^f, 

i^7X« I f'^Mf iOiN NUMBtK: ' f 12//i7X« FACLL FY: aEgulHE^d^fS A.NO COSTS FO 

2K • IN ' , 14/) 

9i)i > UHHAf( iX,2A 1,3 X 'RANK: ' 1 12X.< I' ,9X«2' , 9X« i • , 9X • ' • f^X* SUBT tiF AL • ,6X 

!• TGf ALS' /i2X«FTt^ ' ,'tX,4Fi0.2f l2X,Fi2.2/ 12X • SmLAK 1 1 S « f ^F 1 0 •O ,Fi2,0/ 

.::i2X« SUf^PORT* f'»lX,?,Fi2«0/ ) 
960 FORMAT IVA^^ICOS TABLE: 15- • t 1 4, ^UX • CAT £ : •,^A'r/^ ^U.>l NUMOER: 'fA^t, 

i^7X« I fERAF ION NUiUMiK" • , 1 2 1 / / I OX ' SPAC b: R tUvJ IR bMt£M ( 5 ANU CUSFS FOR 

2 • ,A4,« IM • lU/) 
9oi FUR.SA r(/ lX2A4t I 7X, •CLASSKMS CLASSLA*} KtSEAKCH OFPICE SUPPOH 

IT TUTAi* / UX IMVENTURV* , 3X , 6 i 10/ HX • W ASF • , 3X,ol 10/ IIX* P;<UJ 

^tCT COSrS' ;6: lO.O) 
97J FURMAT ( ' IVA* ICU3 TAtiLE: C«,5JX«0ATE; •,2A4/« HUN ;4U;-<afcR: •,A4,*»9X 

i'YhAHS Irtl^ RUiW^ * » i2//30X • SUMHARY Or VAHlAdLc CUi F -i • // 9X» VE A«< ♦ , 

a« SAl-ARIKS/SU^^POKT' ,-3X'M^U-PLANrSLiX'0t:llT S-^WiCE- 9X«T0FAL»/) 
971 FaK-iAT (9X, I't , 3(^X• ImAI), 5X,14AL) 
9/^ FGKMATi SX, 14 ,3^.Xf I'iAl ,5Xf i^iAl) 

bUeaCUTINL-: FCLEAR(A,rn 
UlKcNSICN A( 1) 
UC ICO 1 = 1 ,N 
lOU A( I) =0* 
RETURN 
lNU 

bLL'KCtriNE NCLEAR(M,N) 
U I^lNSICN Hi I) 
jG iGO l-'lfN 

luo Hi n ^-0 

Kn rUHN 

bUBRCUT INE CGNVI ( N VKR , NO I G, KVAL I 



C NARR IS INrEGER<'2 ARRAY RETURNtJ AS < NO 1 G U NU I 1 ) / 3 ) CHARACiERS 

C CF RES'JLT RtADlNG LEFT TO RIGHT IN AjCcNJiNG J'J.^il^CK IPTb . 

C NUIG IS THt LARGEST NUHbER UF UIGIFS (li4CLU01N,i MiriUS SIGN) THAT 

C MIST EE ASSUMED FOR NVAL. NUiO MUSF bii ^ < N3I0 < 12, 

C KVAL OR FVAL ARE THE NUMtiER TO afc FKANSLArtO. FHlY COME IN 

C THROUGH DIFFERENT ENTRY POINTS 

C NVAL IS ThE VALUE 10 TRANSLAFbOt CHANGED FROM EITHER KVAL 
C CK FVAL 

C NJIGf KVALf FVALf ARE NOT ALTERED dY SGtiRQUriNE* 

C ZiiRU RETUKNS AS tiLANK* ^ 

C CVEHFLOw (EITHER POSITIVE Oa.NEGATIVL) RETURNS AS JN THE 



RlGHTMUSr PQSITJON. 
A ^i.NtS SION I'S KlGHf ADJUSreu rO NUMBcrt, 

r.j LhHOf^ ckeckIng is oone, 

UATA oUN< /'O* r t '^^ »» 3* ^' » •S* . f • 7' »« ^» # • 7* ,«-• t • •f't't 

UATA r. reST /IC llOO, iOOOiiOOOOflOOOQOt lOCOoOO f 100 COOuO, 10 0000000* 

NVAl=KVAL 

'< = Jl)NK( 12) 

iF ir.VAL .LT. 0) K>= JUNKllU 

N-1 . ^> - 

Nl2Nu2 *f»nNJ+l 

UC t j = i,NbNO 
1 WAPKC J>-JUNKil2) 

l^kH'i^L •OE. NiTl5ST(r;L)IGJ ) GU TO 6 

IF (r.VAL .Lt. -NTEST(ND10-U ) GO TO 6 

nVAu -IAJS(NVAI) 
^ uC i O^LfM 

jGKt^fiVAt /N 

IF t JGKC • LU, 0> GO TO 4 

XUKc = (.%VAL /N)*10 
JOKl - JUKd-KOKE 
J WAKKlNE.NOa-J)=JUNK( JOKE+1) 
1F(H .tvJ. NtiNO) GO TO 5 

IF i¥ .or, NENC) M=NENO 
NARK INENL)2-L+1)=K 
IFi^CKE ^EG* 0) GO TO 5 
.^AHK { ^EN02-t + l * = JUNK(13) 

uG rc 2 

6 NARK tNENO) =JUhKi 14) 
v»G fU 5 

tNTrtV CUNVfi(NAf<R,MUlG,FVAL) 
wVAL^FVAL + SlGIU .5f FVAL ) 
olj i L 7 
tNO 



